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We derive a factorization theorem for the Higgs boson transverse momentum (pr) 
and rapidity (Y) distributions at hadron colliders, using the Soft Collinear Effective 
Theory (SCET), for > > A-qcd where denotes the Higgs mass. In ad- 
dition to the factorization of the various scales involved, the perturbative physics at 
the pt scale is further factorized into two collinear impact-parameter Beam Func- 
tions (iBFs) and an inverse Soft Function (iSF). These newly defined functions are 
of a universal nature for the study of differential distributions at hadron colliders. 
The additional factorization of the py-scale physics simplifies the implementation of 
higher order radiative corrections in a s (pr)- We derive formulas for factorization in 
both momentum and impact parameter space and discuss the relationship between 
them. Large logarithms of the relevant scales in the problem are summed using the 
renormalization group equations of the effective theories. Power corrections to the 
factorization theorem in px/n^h an d ^-qcd/pt can be systematically derived. We 
perform multiple consistency checks on our factorization theorem including a com- 
parison with known fixed order QCD results. We compare the SCET factorization 
theorem with the Collins-Soper-Sterman approach to \ow-pt resummation. 
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INTRODUCTION 



The Higgs boson is the last undiscovered particle of the Standard Model (SM), and 
its discovery is a major goal of both the Tevatron and the Large Hadron Collider (LHC) 
physics programs. If a scalar particle is discovered at either collider, the measurement of 
its properties will be crucial to determine whether the particle found is the Standard Model 
Higgs boson, or whether it hints at physics beyond the SM. The theoretical community has 
devoted significant effort to understanding precisely the production cross section and decay 
widths of the SM Higgs particle in order to facilitate such studies, as reviewed recently in 
Ref. pp. The dominant production mode at both the Tevatron and the LHC is the partonic 
mechanism gg — > H proceeding through a top-quark loop [2H6]. The perturbative QCD 
corrections are known through next-to- leading order in full QCD [TIE], while the corrections 
in the m t — > oo limit are known through next-to-leading order |9] and next-to-next-to-leading 
order [101412] . Resummation of logarithmically-enhanced threshold corrections to the cross 
section has been studied [TBTUB] . The inclusion of such theoretical calculations is crucial for 
experimental searches for the Higgs boson, as they increase the predicted cross section in 
the SM by a factor of two at the LHC and by more than a factor of three at the Tevatron. 

The study of differential distributions of the Higgs boson is also needed in experimental 
analyses. For example, for a SM Higgs in the mass range 130 GeV < < 160 GeV, 
one of the most promising discovery modes is through the partonic process gg — )■ h — > 
W + W~ — > i+vi'v. Since the final state contains two neutrinos, reconstruction of the Higgs 
mass peak is not possible. An understanding of the kinematic distributions for both signal 
and backgrounds is needed in this search channel. The NNLO differential distributions in 
the m t — > oo effective theory were obtained in Refs. [TTH2T)] . and detailed studies of the 
effects of experimental cuts on Higgs boson cross sections have been performed [2TH23] . 
However, a large reducible background comes from pp — > ti — > bW + bW~ — > i + vl~v + 
jets. Such backgrounds are brought under control with a series of cuts which include a 
jet veto so that any process involving a jet with high transverse momentum, taken to be 
roughly p T > 20 GeV [23H26] , is rejected. Such a cut selects Higgs boson with primarily low 
transverse momentum, and therefore a proper implementation of such jet vetoes requires 
a good understanding of the Higgs differential distributions at low p? where resummation 
of large pr/m^ logarithms is necessary. The study of the low-transverse momentum region 
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in hadronic collisions has been under investigation since the early days of QCD [2TH29] . It 
has been studied for the Higgs boson following the seminal analysis of Collins, Soper, and 
Sterman (CSS) [301 EI] in several works [32H36] . 

The purpose of this paper is to derive a factorization theorem for the Higgs transverse 
momentum p T and rapidity Y distribution, in the region Q ~ m h 3> p T ^> Aq C d, using the 
Soft Collinear Effective Theory (SCET) [37H39] . Here Q and rrih denote the partonic center 
of mass energy and the Higgs mass respectively. Although we focus on Higgs production, 
our methods and results can be immediately generalized to the differential distributions of 
any one or more color neutral particles. The factorization theorem we derive has the form 

d 2 a 7r 2 f 1 dx\ f 1 dx 2 f 1 dx[ f 1 dx' 2 



dp^dY 4(iV c 2 — 1) 2 Q 2 J Xi J x 2 J Xl x[ J X2 x' 2 

x H(x u x 2 , /iq; {it)G v (x!, x[, x 2 , x' 2 ,p T , Y, /j tT )f i/P (x' l , pjr)fj/p(x' 2 , /xr), (1) 

where Q is the hadronic center of mass energy, H is the hard Wilson coefficient arising 
matching QCD onto SCET, and fi/p is the standard parton distribution function (PDF) for 
taking a parton of species i from the proton. Q 1 ^ is a perturbative coefficient at the pt scale 
that has the form 

G'^x^x'^x^x'^pt.Y^t) = J dt+ J dt n J Jo(b± Pt) gaa9p„ 

X ^n;!^^'^'^'^)^,^-^,^,^,^) (2) 

x l x 2 



x S \x 1 Q -e Y ^pl + m 2 h - -^,x 2 Q -e Y ^Jp 2 T + m 2 - b±, // T ) 

which is a convolution over the collinear functions l^n g % an d the Inverse Soft Function (iSF) 
S^ 1 . Logarithms of rrih/pT are summed by the Renormalization Group (RG) equations in 
SCET and are encoded in H(x\, x 2 , \iq\ ht) which is the hard coefficient evolved from the 
renormalization scale hq ~ m h down to /iy ~ px- The logarithms of Aq C d/pt ar e summed 
via the standard DGLAP evolution of the PDFs and are encoded in the PDFs evaluated at 
Ht ~ Pt- The factorization formula in Eq.([T]) is derived by matching QCD onto a sequence 
of effective field theories EFT: 

QCD(n 7 = 6) ->■ QCD(n/ = 5) ->■ SCET PT ->■ SCET Aqco , (3) 

which is shown graphically in Fig. [l| The first step Q,CD(nj = 6) — > QCD(nj = 5) denotes 
the usual procedure of integrating out the top quark to get an effective coupling of the Higgs 
to gluons. The Higgs production mechanism then proceeds via this effective coupling. The 
hard scale Q ~ is then integrated out by matching onto SCET PT , which describes the 
dynamics of soft and collinear modes with transverse momenta of order px- The factorization 
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theorem in SCET PT takes the form 
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dp 2 T dY ~ 4(iV c 2 -l) 2 J dp * dp ~ h j ^ J (2^ 2 



d 2 b 



X 5 



P h ~ eY V 'Pt + m l 



5 



Ph- e ¥ \/PT + m h 



X 



ptph - k h± - m h 

dxi J dx 2 J dt+ J dt^H(x!X 2 Q 2 , i^q] B^ix^t^^x, fx T ) B naP (x 2 ,t^,b±, fi T ) 



t- t + 
x S {xiQ-p h --^,x 2 Q-p'l--^,b ± ,fi T ), 



where the collinear functions B^ n are the Impact-parameter Beam Functions(iBFs). The 
iBFs B^f n are extensions of the beam functions that appear in [^CTf BT] and reduce to them for 
b± = after contraction of the transverse indices a and (3. The beam functions of Ref. [H] 
were shown to have wide applicability to the analysis of observables at the LHC. The iBFs 
are proton matrix elements evaluated at the scale \xt ~ pr- The iBFs are matched onto the 
standard QCD PDFs by performing an OPE in Aqcd/pt and the logarithms of Aqcd/pt 
are summed via the standard DGLAP equations used to evaluate the PDFs at the scale 
fiT ~ Pt- This is shown schematically in Fig. [T]and gives the final form of the factorization 
theorem shown in Eqs.([l]) and (2) where the collinear functions T^ n . g i are just the iBF to 
PDF matching coefficients. 

While the factorization and resummation of transverse-momentum distributions has been 
studied extensively in the QCD literature, and SCET analyses (121 H3] have been performed 
in the past, our analysis contains several interesting differences that we believe are worth 
further investigating. A summary of the main points of this paper is given below: 

1. We derive a factorization theorem for the Higgs transverse momentum and rapidity 
distributions using effective field theory methods. A clear separation of the dynamics 
associated with the scales Q ~ m h ^> p T ^> Aq C d into perturbative Wilson coefficients 
and standard QCD PDFs is achieved. Large logarithms of ratios of the relevant scales 
are summed using RG equations in the effective theories. Power corrections in pr/n^h 
and Aqcd/pt can be systematically derived by going to higher orders in the power 
counting of the effective theories. 

2. In addition to the factorization of the scales rrth > pt > Aqcd, the perturbative 
physics of the px scale is further factorized into an iSF S^ 1 and two distinct collinear 
functions 1^. gi and 1^ gi - This additional factorization simplifies the structure of higher 
order perturbative corrections at the pt scale. They can now be obtained through 

raP n „A T a P 



higher order computations of the simpler perturbative functions S 



The factorization in SCET naturally occurs in terms of purely collinear PDFs and 
soft functions. The purely collinear PDFs differ from the standard QCD PDFs by soft 



(4) 
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FIG. 1: Structure of the factorization theorem. The QCD (nj = 5) theory obtained after integrat- 
ing out the top is matched onto SCET PT at the fi ~ scale followed by RG running down to 
the [i ~ pt scale summing logarithms in rrih/pT in the process. Using the soft-collinear decoupling 
property of the leading order SCET PT Lagrangian, the cross-section is factorized into a n-collinear 
iBF (Impact-parameter Beam Function), a n-collinear iBF, and an iSF (Inverse Soft Function). The 
iBFs are then matched onto the standard QCD PDFs at the fi ~ pt scale and the logarithms of 
Aqcd/pt are summed via the DGLAP equations which determine the PDFs at the [i ~ pt scale. 



subtraction terms required to avoid double counting soft emissions that are already 
contained in the soft function. However, the equivalence of soft zero-bin and soft 
subtractions allow us to rewrite the factorization in terms of the standard QCD PDFs. 

We give expressions for the factorization formula in both impact-parameter space (b- 
space) and momentum space and discuss how they are related. The factorization can 
be formulated entirely in momentum space; see Eqs.Q, ( [53| , and (52). However, the 



matching coefficients X".^ and l^f ■ are obtained by matching the iBFs in b-space onto 
the standard QCD PDFs. 

5. The iBFs B^ n {x-, t, b±, /i) that appear in the SCET PT factorization theorem in Eq.(4) 
are extensions of the beam functions introduced recently in |10j E] . The additional 
functional dependence on b± found here is required to facilitate resummation of the 
low pt region, and our functions reduce to those studied previously for b± = 0. 

6. The iBFs are similar to transverse- momentum dependent PDFs studied previously in 
the literature [HH17] . except for an additional dependence on a residual light-cone 
momentum component. The iBFs are defined in covariant or non-singular gauges 
and are invariant under non-singular gauge transformations. In singular gauges, like 
the light cone gauge where the gauge potential is non-vanishing at infinity, the iBFs 
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must be modified [HHU] with a transverse gauge link to be gauge invariant under the 
singular gauge transformations. However, the convolutions over perturbative quan- 
tities ® B na p ® S" 1 is fully gauge invariant since it is obtained by applying a 
soft-collinear decoupling transformation to a manifestly gauge invariant operator in 
SCET PT . Correspondingly, after matching the iBFs onto gauge invariant PDFs, the 
quantity Q % i = 1^\g,i ® <8> is also gauge invariant. As a result, even though 

the iBFs are defined in covariant gauges, the factorization theorem is completely gauge 
invariant as expected. The additional functional dependence on the residual light-cone 
momentum avoids singularities that occur when this variable is inclusively integrated 
over, as it is when considering the usual transverse-momentum dependent PDFs. 

7. The iBFs come with perpendicular indices a, (3 which are contracted between the n- 
collinear and the n-collinear iBFs. This structure has not been discussed in previous 
SCET analyses. This non-trivial index structure encodes dot products between the 
momenta of final state gluons emitted in different directions, and is needed to sepa- 
rately calculate contributions from n and n gluons. 

8. The Landau poles associated with impact-parameter integrations in the CSS formu- 
lation can be avoided in the effective field theory approach employed here. As al- 
ready mentioned, the factorization theorem can be formulated entirely in momentum 
space with no reference to any impact-parameter or impact-parameter integrations. 
The perturbative function Q % i and the PDFs in Eq.Q are naturally evaluated at the 
renormalization scale /i ~ Pt so that one never encounters the Landau pole associated 
with the limit \i — > 0. 

9. The iBFs are not specific to Higgs production, but have universal applicability to the 
study of differential distributions in gluon-initiated processes at colliders. Analogous 
functions exist for quark-initiated processes, and can be derived following the approach 
discussed in this work. 

10. In the non-perturbative region px ~ ^-qcd, where the effective field theory power 
counting in Aqcd/pt breaks down, a phenomenological model can be introduced for 
the quantity C/ u analogous to what is done in the standard approach. 

11. We perform several consistency checks on the factorization theorem including a com- 
parison with known fixed order QCD results. We ensure that the renormalization scale 
dependence properly cancels among the various factored objects in the factorization 
formula, which is a non-trivial check on the structure of the factorization theorem. 

There have also been many other recent [HJ 0HH53J I53H57] developments in the application 
of SCET to describe other processes at hadron colliders and it is a promising avenue to 
pursue in the LHC era. 
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Our paper is organized is follows. In Section [II] we review the calculation of Higgs produc- 
tion through the gluon-fusion mechanism in perturbative QCD. We derive our factorization 
formula in Section III We discuss the factorization formula and compare to other results in 
the literature in Section IV In Section [V] we calculate the various quantities in fixed-order 
perturbation theory, while Section VI is devoted to resummation of logarithms. We per- 
form a series of consistency checks on our factorization formula including a comparison with 



fixed-order perturbative QCD in Section VII Finally, we conclude in Section VIII 



II. HIGGS PRODUCTION IN QCD 



We begin by reviewing the gluon-initiated production of a Higgs boson in QCD. The 
coupling of the Higgs boson to gluons arises primarily from a top-quark loop. For rrih < 2rn t , 
we can integrate out the top quark to derive an effective coupling of the Higgs boson to 
gluons [21 El El EE] • The effective Lagrangian is given by 

£ mt = C GG h - G^ v G% v , C GGh = jl + — + 0(a^)j , (5) 

where C GG h, the Wilson coefficient in the MS scheme, is known through 0(a 4 s ) [T31 1591462] . 
Calculations of the total cross section at higher orders in QCD perturbation theory using this 
effective Lagrangian have been shown to reproduce the result of the full theory to percent- 
level accuracy when rrih < 2m t if the effective-theory cross section is normalized by the full 
top-quark mass-dependent leading-order cross section [8, 63J. Additional corrections to the 
transverse momentum spectrum of 0(p^/mf) are also present in the full theory. These affect 
the low-pr region at the percent level, and have recently been studied in Refs. [64, 65J. An 
effective field theory approach for studying such corrections at higher orders has recently 
been developed [66, 67]. The scale \i at which the coupling constant is evaluated should be 
chosen as fi ~ m t to minimize logarithms that appear in the 0(a^) expression for C GG h- 
For notational ease we define 

g 2 c=-4C GGh (6) 

and express later results using c. 

The Higgs boson must recoil against at least one parton in order to have non-zero trans- 
verse momentum. At leading-order in perturbative QCD, three partonic processes contribute 
to Higgs production at non-zero pr- gg —> hg, q(q)g —> hq(q), and qq —> hg. We focus here 
on the dominant process gg — > hg. The other partonic production channels are suppressed by 
smaller PDF luminosities, and can be straightforwardly included in the analysis if desired. 
The diagrams contributing to this process are shown in Fig. [2] We denote the incoming 
proton momenta as pi,P2, the outgoing parton momentum as p^, and the outgoing Higgs 
momentum as ph- We define the light-cone vectors n and n through 

n» = (1,0,0,1), n" = (1,0, 0,-1). (7) 
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FIG. 2: The leading order diagrams in QCD that contribute to Higgs production with non-zero 
PT- 

The momenta of the initial hadrons and the Higgs boson are given by 



where Y — 1/2 ln(n • pt/n ■ ph) is the rapidity of the Higgs boson. The partons entering the 
hard-scattering process have momenta and virtualities given by 



where are the usual Bjorken momentum fractions. 

Several distinct kinematic regions contribute to production of a Higgs boson at low trans- 
verse momentum. In the first, the Higgs recoils against a n-collinear gluon with the following 
momentum in light-cone coordinates: 



This corresponds to a Higgs boson produced at high rapidity and low transverse momentum. 
A similar region exists with the Higgs recoiling against an n-collinear gluon with scaling 
P3 ~ fnh(^,PT/ m hyPT/ m h)- Finally, the Higgs boson may recoil against a gluon with an 
"soft" momentum p 3 ~ (pt,Pt,Pt)- This corresponds to the production of a Higgs at central 
rapidity and low transverse momentum. Since the final-state gluons in the low-momentum 
region are restricted to be either soft or collinear to the initial partons, large logarithms 
of the form In (m^/pr) appear in the perturbative expansion and must be resummed to 
all orders. The n-collinear, n-collinear and soft gluon modes will be used to construct the 
effective theory that facilitates this resummation. 

For future use we reproduce here the differential and total cross sections for the gluon- 
gluon scattering process. The hadronic cross section can be expressed as a convolution of 
parton distribution functions and partonic cross sections: 





Pi = XlPl + O(Aqcd), P2 = X 2 P2 + O(Aqcd) 






(10) 




(11) 
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where s, t, and u are the usual partonic Mandelstam variables. For production of the Higgs 
with non-zero Pt, the differential partonic cross section is given by [68] 

da 7i / a s \ 3 J m\ + s A + t 4 + u A 
~dl ~ 384w 2 WJ { §fa 

The total partonic cross section for gg — >• h through next-to-leading order in QCD pertur- 
bation theory is E] 




where z = m 2 h /s. This result assumes the scale choice fi 2 = s. The dependence of the 
partonic cross section on the renormalization and factorization scales can be restored by 
using the known renormalization group running of the cross section. The result is presented 
in Ref. [9]. 

III. EFT FRAMEWORK 

We derive a factorization theorem via a sequence of effective theories 

QCD(ra/ = 6) ->■ QCD(n / = 5) ->■ SCET PT ->■ SCET Aqcd , (14) 

which factorize the physics associated with the different scales Q ~ rrih ^> Pt ^qcd into 
calculable perturbative functions and standard QCD PDFs. As we are assuming that the 
mass of the Higgs is sufficiently small {m h < 2m t ), we can integrate out the top quark in the 
matching step QCD(n/ = 6) — > QCD(nj = 5) to obtain an effective coupling of the Higgs 
boson to gluons. The cross sections obtained using this effective theory were described in 
Sec. [TTJ To derive a renormalization group equation allowing resummation of large logarithms 
In (nth/pr) that appear at low transverse momenta, the matching to SCET PT is required. 
The soft-collinear decoupling property of the leading order SCET PT Lagrangian also leads 
to a factorization of the soft and collinear sectors, which simplifies calculations of the cross 
section in the low px region. Finally, the matching to SCETa qcd expresses the cross section 
in terms of the standard parton distribution functions. We describe in this section the details 
of each stage in the matching in QCD(nj = 5) — > SCET PT — > SCETa qci? . 

A. QCD to SCET PT 

As already mentioned, the perturbative expansion in QCD for the transverse momentum 
spectrum of the Higgs contains logarithms of rrih/pT. In the low transverse momentum 
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region Aqcd <C Pr ^ rnh, these logarithms become large and must be resummed to all 
orders in perturbation theory. In the effective theory formulation, this is done by matching 
QCD onto the effective theory SCET PT , which describes the dynamics of the degrees of 
freedom recoiling against the Higgs, and solving the RG equations of the effective theory 
operators. The effective theory SCET PT is formulated in terms of collinear and soft modes 
with momentum scalings 

Pn ~ m h (ri 2 ,l,Ti), p n ~ m h {\, r] 2 ,r]), p s ~ 77^(77,77,77), rj ~ — , 

(15) 

where p n ,Pn, and p s denote typical momenta for the n-collinear, n-collinear and soft modes 
respectively. The effective theory has a well defined power counting in the parameter 77 and 
has distinct quark and gluon fields for each of these modes. We note that ultrasoft gluons 
with scaling p us ~ mh(i] 2 , i] 2 ,r] 2 ) would not produce a Higgs boson with sufficient transverse 
momentum, and therefore do not need to be included as separate modes in the effective 
theory. The gluon fields A^p n (x), A^p_(x), and A^Jx) destroy n-collinear, n-collinear, and 
soft gluons respectively. The presence of distinct collinear and soft gluons requires the theory 
to be invariant under collinear and soft gauge transformations EH] ■ The momenta of the 
effective theory modes are separated into label p and residual k parts 

p**=pP + fcM p^ ~ rn h (l,rj), ~ m h 7] 2 . (16) 

Derivative operators are similarly separated into label and residual operators so that, for 
example, a derivative acting on the n-collinear field takes the form 

fit 1 _ 

^4yP + Pj+i9", (17) 

such that the label operators act on the label momentum subscripts 

VnKfiM = n ■ pA^ n (x), = PX<p», (18) 

and the residual derivative operator acts on the residual co-ordinate dependence x M . We note 
that such a field with label momenta can be written explicitly as a Fourier transform of a 
standard quantum field. As an example, a field with no dependence on residual coordinates 
can be expressed as 

X f M = J ^e-^ 2 X(y). (19) 

As already discussed, after integrating out the top quark, the gg — > h process is mediated 
by the effective QCD operator 

r i An 

Oqcd = g 2 hTr ]G^Gr\ = (20) 
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where C mt is given in Eq. In SCET PT , the leading order operator that mediates this 
process is (T5J H2] 



0(u l ,u 2 )=g llv hT{Tr S^gB^J^SlS^gB^^Sl } 



(21) 



where we have suppressed the perpendicular labels on the soft and collinear fields, the sum 
of which are constrained to be the negative of the Higgs transverse momentum, and we use 
the standard notation 



(X) u = XS(tu - V 



(22) 



and the big brackets in Eq. (21) denote the fact that any label operators appearing inside 



do not act outside the brackets. The B field is defined as [70 

_ ' n 

and the soft Wilson line S n in position space is given by 

S n ( x ^) = Pexp (ig s [ dsn- A a s (x» + sn») 



(23) 



(24) 



with an analogous expression for S n . The covariant derivatives are dressed with momentum 
space Wilson lines so that 



WJx) 



E Ex p( 



perms 



p-n ■ A n>q ) 



(25) 



and are the usual covariant derivatives 

in ■ D n = P n + gn ■ A UjP , iD^ = V» x + gA^~, in ■ D n = in ■ d + gn ■ A n}P . (26) 
We match the QCD operator onto the the effective SCET PT operator via 



Oqcd = / dui I du 2 C{ui,u 2 )0{ui,u 2 ) 



(27) 



and determine the Wilson coefficient by computing the gg — )■ h process on the LHS in QCD 
and on the right side in SCET PT . At tree level the Wilson coefficient is 



C (( Vi,u; 2 ) 



C U\UJ2 



V 



(28) 



The process gg — > gh where the transverse momentum of the Higgs is of order mhf] ~ Pt is 
mediated through an extra soft or collinear emission from the Wilson lines in 0(u>i, u 2 )- 

We note that one can instead implement a two-step matching procedure by first matching 
onto an intermediate SCET with scaling parameter £ ~ \jvTl m h- The relevant modes in this 
theory are collinear and ultrasoft gluons with momenta p n ~ m h (£ 2 , 1, £), p n ~ m^,(l, £ 2 , £), 
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and p us 
O ~ Tr 



Y B 

1 n J -'n. 



ytyi B- 1 fJ 



The leading operator in this theory takes the schematic form 
, where Y n> n are ultrasoft Wilson lines built up from eikonal 
interactions with the collinear gluons. The ultrasoft gluons of the intermediate theory just 
become the soft modes of SCET PT , and we replace Y nfi — > S n)fi to obtain Eq. (21). The 
Wilson coefficient from matching these operators is unity to all orders as also seen in |42j. 
A more detailed explanation of this point can also be found in [71] . 

For convenience, in deriving the factorization theorem, we will work with soft fields that 
are in position space and collinear fields that are in position space conjugate to the transverse 
label momentum so that no label momenta of the transverse momentum appear in the soft 
and collinear fields. This amounts to not separating from momentum components of order 
myj] a residual part of order m^ff '■ The separation into label and residual components is 
employed only for the hard light cone collinear momentum components of order m^. 



B. Factorization in SCET PT : iBFs and the soft function 

It is easier to work in terms of the hadronic Mandelstam variables u and t instead of 
Pt and Y, which correspond to the Higgs transverse momentum and rapidity respectively. 
These two sets of variables are related by 



u = (p 2 - Ph) 2 = rn 2 - Q Jp T + m? e Y , 

t = (pi ~ Phf = m 2 -Q sJp 2 T + m 2 e~ Y . (29) 

The transformation between these sets of variables has a rather simple Jacobian given by 

dudt = Q 2 dp 2 T dY. (30) 

Thus, a restriction on the u and t Mandelstam variables is equivalent to a restriction on the 
Pt and Y of the Higgs. The double differential cross-section in the Mandelstam variables 
can be written in SCET as 



d 2 a 1 



1 



d 2 p h± f dn-p h dn-p h _ _ 2 

(2n)9(n-p h + n-p h )b{n-p h n-p h -p h -m h ) 



dudt 2QH4J J (2vr) 2 J 2(2tt) 2 

x 5{u- (p 2 - Ph) 2 )5(t - (p a -p h ) 2 ) Yl J2\ C ^ U2 ^ (hX n X n X s \0(u u u 2 )\pp) | 2 

initial pols. X 

x [2^\ Pl +p 2 - P Xn - P Xn - P Xs - p h ), 

(31) 

where O and C denote the SCET PT operator and the matching coefficient respectively. The 



® symbol denotes a convolution in the label momenta u% t 2 as in Eq. (27). Note that the 
constraint delta functions 8{u — (p 2 — Ph) 2 ) an d 5(t — (pi — Ph) 2 ) restrict the final states 
to those that satisfy u = (p 2 — Ph) 2 and t — (pi — Ph) 2 , or equivalently pick out the states 
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with the corresponding values of pr and Y. The states X n , X n , X s correspond to final state 
particles with the n-collinear, n-collinear and soft momentum scaling respectively. It is only 
the states with such momentum scalings that will have a non-zero overlap with the SCET PT 
operator 0(ui,u 2 ). The overall factor of 1/4 in square brackets in Eq. (31) is from the 
average over the initial proton spins. 

Using the fact that the soft and collinear modes are decoupled in the leading order 
SCETp^ La grangian we arrive at the factorization formula 



d 2 a 



(27T) 



du dt (N 2 - 1) 2 8Q 2 



dn ■ p h / dn ■ p h / d 2 k^ / dk^d 2 k^ / dk n d 2 k^ / d 4 k 



dx d 2 x± f dy d 2 y± f d 4 z 



2(2vr) 



2(2tt) 



(2tt)< 



\k+x -ik^-x ± c ^k n y+-ik±-y 1 _ ik s -z 



x 5 (u - m\ + Qn ■ p h ) 5 (t - m\ + Qn ■ p h ) 5 (n ■ p h n ■ p h - k 2 ± - m 2 ^ 
x J dutd^lC^,^, i^)\ 2 J^(ui,x~ ,x ± , /j) J na ^(u 2 ,y + ,y±^) S(z, /j) 
x 5 ( Ul - n - Ph - k^ - k~) 5(u 2 -pt-K- K)5 {2) {ki + k^ + k^ + ki 

where the jet and soft functions are defined as 
jf(u u x-,x ± ,fi) = ^ (pi\[9Bt f ^{x-,x ± )5(r-u 1 )gB^(0)]\pi^ 



(32) 



initial pols. 

Jf(ui u y + ,y ± ^) = (P2\[gBt±p(y + ,y±)5(V - u 2 )gBi ±a (0)]\p 2 ) 

initial pols. 

" T b ^„T c S't^T D 4j (0)J |0). 

(33) 



S(z,fx) 



\T [Tr (s n T D SiS n T c Si) (z) 



T is the time-ordering symbol, and T denotes anti-time ordering. Details of the derivation 
of this formula are given in appendix |A} The above factorization theorem can be brought 
into a more concise form involving a simpler convolution structure so that 



d 2 a 



(27T) 



du dt (N 2 - 1) 2 8Q 2 



dPhd>Ph d k h dio\du 2 



db + db~d 2 b 



4(2tt) 



± e i { u 1 - p -)b+ p ^ 2 -p+)b 



e 2 ' 



x [u-m 2 h + Qp h ] 5 [t - m\ + Qp\\ 5 \p+p h - k 2 h± - m 2 h 

X \C(tOi, U}2,/J t )\ 2 Jn P (^l,b~,b±,p,) Jnal3(u2,b + ,bj_,fl) S \b + , b~ , bj_, (i) . 



(34) 



We recast this factorization theorem in terms of jet and soft functions that have momentum 
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space light cone coordinates as 
<Pa (2tt) /■,+._ /\ 2 ,± /• 



dudt (N* - 1) 2 8Q 2 J r h j n J 

x 5 [u - m\ + 5 [t - m 2 h + Qp+] 5 p^-pl ~ k 2 h± - m\ J du} 1 du2\C(u ll u 2 ,fM)\ i 
x J dk+dk^ Bf(u 1 ,k+,b ± ,fi) B^(uj 2 , fer, 6 ± , /i) -p^ - k n ,u 2 -p+ - A;+,6j 



where we have defined the hybrid-fourier space jet and soft functions as 

Bf(u 2 ,kr,b ± ) = j '^-e^jf^b+^ri, 



( 



S(Co u Co 2 ,b ± ,f,) = I d ^^ e ^ b+ e^ b -S(b+,b-, 6 ±)A i). 



(36) 



The hybrid jet functions B^ n (uj, k ± ,b±, /x) are similar to the functions that appeared in [3D] 
and more recently in [31], but differ because of their dependence on the impact parameter 
b± and because their perpendicular indices a, (3 are not contracted with each other. We will 
refer to these functions as impact-parameter Beam Functions (iBFs) in analogy to the Beam 
Functions of [31] which have b± = 0. These iBFs are implicitly defined with a zero-bin [72] 
subtraction in order to avoid double counting the soft region already present in the soft 
function S(oji,oj 2 , b±,fi). For clarity we will refer to the zero-bin subtracted iBF as a purely 
collinear iBF. We will refer to the iBF defined without a zero-bin subtraction as the naive 
iBF or simply the iBF 1 when the context is clear. These purely collinear iBFs that appear in 
the factorization theorem are in general gauge dependent quantities. This is seen from their 



Note that the iBF is still always implicitly defined with an "ultra-soft" zero-bin subtraction [41] to avoid 
double-counting between the collinear and the ultra-soft modes with scaling p us ~ M(rj 2 ,i] 2 ,r) 2 ). In 
the paper, the phrase "zero-bin subtraction" refers to the "soft" zero-bin subtraction, to avoid double- 
counting between the collinear and soft modes, unless otherwise specified. The ultra-soft and soft zero- 
bin subtractions are distinct since the soft modes are constrained by pt while the ultra-soft modes are 
unaffected. The soft zero-bin subtraction in the beam functions can be recast in terms of the two iBFs 
and an inverse soft function (iSF) as shown in the next section. Now the iBF is defined without the soft 
zero-bin subtraction which is instead taken into account by the presence of the iSF. However, there is still 
a double-counting left over between the collinear and the ultra-soft modes (which are unaffected by pt 
constraints) which is avoided by the ultra-soft zero-bin subtraction in the iBF. This ultra-soft zero-bin 
subtraction in the iBF should be understood to be implicitly part of its definition. Thus, at b± = 0, 
the iBF reduces exactly to the beam functions in |41j . In pure dimensional regularization, the ultrasoft 
zero-bin subtraction is scaleless and vanishes. We thank I.W. Stewart and F. Tackmann for detailed 
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dependence on the impact parameter b± ^ which leads to a spatial separation between 
the fields in the matrix element not connected by a Wilson line. However, this additional 
gauge link can be placed at infinity along the light-cone, and it does not contribute in 
covariant gauges where the gauge potential vanishes at infinity. The iBF is thus well-defined 
in covariant gauges. This is similar to what occurs for transverse-momentum dependent 
PDFs in QCD [HJ-HT]. In light-cone gauge, this additional gauge link at infinity is required 
due to the asymptotic behavior of the gauge potential. It is possible that Glauber modes 
are responsible for building up this extra contribution in SCET [17]. We note that the 
total convolution over the hard Wilson coefficient \C{u\, u 2 , /i)| 2 , the purely collinear iBFs, 
and the soft function in Eq. (35) is just equal to the total perturbative cross-section for 
gluon-initiated Higgs + multi-parton production and thus gauge independent as required. 



C. Equivalence of zero-bin and soft subtractions 

The purely collinear iBFs B^n(u>, k^, b±, fi) defined with a zero-bin subtraction can be 
written as 

K P n(u,k ± ,b±,iJ,) = I3^ n (u,k ± } b ±} fi) - B^ om (u,k ± ,b±,fi) (37) 

where Bnfi( UJ ^k ± ,b±, /1) is the naive iBF or simply the iBF defined without a soft zero- 
bin subtraction but still with an implicit ultrasoft zero-bin subtraction. The functions 
B{no m}( u > k- 1 ' Z 1 ) denote the soft zero-bin limit of the iBFs. It has been demonstrated 
in several processes [73H7H] that the soft zero-bin subtraction is equivalent to a subtraction 
of the soft function. The same holds true in this case allowing us to recast the factorization 
in terms of the iBFs as opposed to the purely collinear iBFs. In appendix [B] we show that 
the convolution over the purely collinear iBFs and the soft function can be replaced with a 
convolution over the naive iBFs with an Inverse Soft Function (iSF) so that 

J du>idu) 2 \C(u)i,u}2, yu)| 2 J dk+dk^B% ts (u)i,k+,b±, /i) B Ra p(u 2 ,k^,b±, /i) 
x S(ui -vl - K,u 2 -pi - k+,b±,fi) 

= J du 1 du 2 \C(u 1 ,u 2 ,fi)\ 2 J dk+dk^B^(u 1 ,k+,b ± ,fi)B nal3 (u 2 ,k^,b ± ,Li) 
x S~Vi ~Ph ~ k n^2 -pt - k+,b±,fi). 

(38) 

It is useful to recast the factorization theorem in the latter form in terms of iBFs and an 
iSF and eventually the iBFs will be matched onto the standard QCD PDFs via an operator 



discussions which led to an understanding of this point while at the Aspen Center for Physics summer 
workshop on Forefront QCD and LHC Discoveries, 2010. We also thank A. Jain, M. Procura, and W. 
Waalewijn for an inquiry which prompted us to explain this point here explicitly. 
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product expansion. The factorization theorem in terms of the iBFs and an inverse soft 
function takes the form 



d 2 a 



(27T) 



dp+dp h / d k h 



d 2 bi 



e -ik£-b ± 



dudt 8Q 2 (N 2 - l) 2 J ^ m j 11 J (2n) 2 

x 5 [u - m\ + Qp~] 5[t-m 2 h + Qp+] 5 p\pl - k 2 h± - m\ J du 1 du 2 \C( y u 1 ,u 2 , fi)\ 2 

x J dk^dkl Bf(u 1 ,k^b±,fx) B na p{u 2 , A£, b±, S~ 1 {u 1 -p h - k n ,u 2 -p\ - k+,b±,fi). 

(39) 



For later convenience we switch to the variables 



Xi = 



n-pi Q n-p 2 Q 

so that the factorization theorem takes the form 



u 2 



00 2 



t n — Qk n , t n — Qk n , 



(40) 



d 2 a 



(27T) 



dp+dp h / d k h 



d 2 b. 



dudt 8(iV c 2 -l) 2 7 * h * h ) h J (2tt) 2 

x S [u - m 2 h + Qp~] 5[t-m 2 h + Qp+] 5 p+p^ - k 2 h± - m\ 



x 



J dxi J dx 2 J dt+ J dt n H(xix 2 Q 2 ^ Q -^ T ) B^( Xl ,t^,b±, ji T ) B nal3 (x 2 ,t n ,b ± , fx T ) 



tt 



x S 1 (x 1 Q -p h - ^-,x 2 Q - pi - y-,b±,fi T )- 



Q 



Q 



(41) 



We have defined 



B?( Xl ,t+,bj_,ri = B?( Xl X,b±,l*), 
Bf(x 2 ,q,b ± ,n) = Bf(x 2 ,k^,b ± ,»), 
H(xiX 2 Q 2 ,^) = \C( Xi x 2 Q 2 ,ij,)\ 2 , 



(42) 



and H(xix 2 Q 2 , hq; /it) denotes the result of RG evolving the function H(x\x 2 Q 2 , //) from 
the scale hq ~ rrih to the scale n ~ pr- In the above equation, the iBFs are written in terms 
of tnn instead of k^n to make it manifest that the iBFs actually depend on t^n, as demanded 
by reparameterization invariance. The choice of the scale \xt ~ Pt will become manifest once 
we perform the Higgs phase space integrals and rewrite the u and t Mandelstam variables in 
terms of pr and Y. We will do this in the next section. The RG evolved H(xix 2 Q 2 , fiQ] ht) 
hard function sums up logarithms of m^/p^. The iBFs are proton matrix elements and will 
give rise to logarithms of Aqcd/pt in the perturbative cross-section that must be resummed. 
For this reason, as discussed in the next section, the iBFs will be matched onto PDFs and 
the logarithms of Aqcd/pt will be resummed via the standard DGLAP evolution equations. 



18 



D. iBFs to PDFs 



The matching of the iBF onto the PDF is given by 



i=g,g,q 



1 dz' 



(43) 



where X^(^-, t^,b±, /x) is the matching coefficient and the gluon pdf is defined as 



—zn ■ pi 



£ (Pil [Tr{^(0)5(P -zn. pi)S ±M (0)}] 



spins 



so that the leading order perturbative expression is normalized as 



(44) 



(45) 



A matching equation analogous to Eq. (|43l) holds for the n-collinear iBF B^. Note that 



in the iBF to PDF matching in Eq. (43), the iBF can match onto quark PDFs beyond tree 



level. In the initial analysis presented in this paper, we ignore this effect. It is simple and 
straightforward to include the effects of the quark PDFs if desired. By noting the the PDF 
is scaleless and that the infrared structure of the iBF and PDF match, as discussed further 
in Section [V] one obtains the following all orders expression for the Wilson coefficient 

= -z[Bf(^z'4,b ± ^) 

We explicitly check this expression at one order beyond tree level in Section [V} An analogous 
expression holds for the n-collinear Wilson coefficient X^ a gi . 



finite part in dim-rcg 



(46) 



Using Eq. (43) in Eq. (41) we arrive at the factorization theorem 



(Pa 
du dt 



(2tt) 



i(N* - 1) 2 Q 2 



d Ph d Ph / d k 



d 2 b A 

W) 



PhPh 



m. 



x 



8[u-ml + Qp- h ]5[t-ml + Qpt} f ^ f ^ f ^ f dx ^H(x l x 2 Q 2 , /i Q ; n T ) 

J0 x l Jo x 2 Jxi x l Jx 2 X 2 

/ dt n dt n 9aa ffjju X n^g,i ( ~T > *n ' & -L > Vt) Z™ g j {^T,t^,b L ^ T ) 
J x l x 2 



t- t + 

x S \x!Q -p h - -^,x 2 Q -pt - -^,b ± ,^ T )f i/ p( y x' l ,fX T )fj/p( x 2^T) 



Next we perform the Higgs phase space integrals to get 



(47) 



(Pa 
du dt 



7l~ 



dx\ f dx 2 f dx^ f dx<2 



4Q 4 (iV c 2 - l) 2 J xi J x 2 J X1 x\ J X2 x 2 



x 



x H(x 1 x 2 Q 2 ^q]^ t ) G lJ (x 1 ,x' 1 ,x 2 ,x' 2 ,u,t,jj T ) f i/P (x' 1 ,iJ lT )f j /p(x' 2 ,iJ lT ), 



(48) 
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where we have defined the u and t dependent function 

d 2 b> 



g iJ (x 1 ,x' 1 ,x 2 ,x' 2 ,u,t,iJ lT ) 



dtl / dt n 



(2tt) 



(m 2 -u){m 2 -t) 



v , Tpa / x l .+ 7 \ -7-/80 j \ 



.r 



3j 9 



m 2 — u t. 
Q Q 



n ^ 



mi - 1 G 



g g 



b ± ,fi T ). 



(49) 



In terms of the and K variables, related to the Mandelstam u and t variables as in 



Eq. (29), we have 
d 2 a 



7T~ 



dx\ f 1 dx 2 f 1 dx[ f 1 dx'r 



X2 ""2 



dp T dY 4(iV 2 - l) 2 g 2 J xi j x 2 J Xl x\ 

x H(x 1 x 2 Q 2 , hq\ ^ t )Q %3 (xi, x[, x 2 , x' 2 , p T , Y, fj lT )f i/P (x' 1 , ii T )f j/P (x 2 , fx T ), 



(50) 



where we have defined the pt and Y dependent perturbative function 

d 2 bi 



G l3 (xi, x' 1: x 2 , x' 2 ,p T , Y, ht) 



dtz / dt 



n I n 



(27T) 



;M\ b ±\pr) 



x 5- 1 (x 1 g-e y ^/p2 +m 2 - -A,x 2 g-e- y ^/p 2 . + m2 - 



G 



E. Momentum space vs impact-parameter space 

As seen in the last section, the factorization formula for the px, Y and u, t dis- 



tributions involved the functions G l3 {x\, x[, x 2 , x' 2 ,pr, Y, //y), defined in Eq. (51), and 



G %3 (xx, x[, x 2 , x' 2 , u, t, pLr), defined in Eq. (49), respectively. These functions are defined 



in terms of impact-parameter space Wilson coefficients X^ a n -(z,t, b±) and the impact 



,n;g,i 



parameter space iSF S l (k ,k + ,b±). We can now introduce momentum space Wilson 
coefficients and a momentum-space iSF via 



jl 3a (7 th, 



S- x {k-,k + ,k } 



d^b 
Air 2 
d 2 b± 
4^ 



—ik 1 -b 1 c — 1 



s- L (k-,k + ,b ± ), 



(52) 
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so that the function Q^(xi, x[, x 2 , x' 2 ,Pt, Y, A*t) can be written as 



g ij (x 1 ,x' 1 ,x 2 ,x' 2 ,p T ,Y, l , T ) ~- - I <!tt I dt- I s fk : r I <Ht I <l 2 k^ 5ij>T \ k n+ k n+H\) 



Pt 



1 

T l3a i x \ ,+ il \ q-fla / x 2 ,_ , _|_ \ 

(53) 



with no reference to The expression for G l -*{xi, x[, x 2 , x 2 , u, t, /it) can be trivially obtained 



from Eq. (53) by rewriting the px and Y variables in terms of the u and t variables by 



using Eq. (29). Eq. (50) together with Eq. (53) gives the factorization theorem entirely in 



momentum space with no reference to the impact parameter. This form of the factorization 
makes manifest the px dependence and how it is related to the momenta of the final state 
particles involved. 

However, the impact-parameter space formulation is necessary for the matching of the 



iBF onto the PDF. This matching must be done in impact-parameter space as in Eq. (43). 



It is only in impact-parameter space that the infrared singular structure of the iBFs and 
the PDFs becomes manifest. Thus, in order to obtain the momentum space factorization, 
one must first go to the impact-parameter space formulation and obtain the iBF to PDF 
matching coefficients I^ n . g i (z,t,b±) and then fourier transform to momentum space. 



IV. FACTORIZATION THEOREM: SUMMARY AND DISCUSSION 



A. General discussion 



In this section we summarize the main points and give a more detailed discussion about 
the structure of the factorization theorem. We also compare our results with both previous 
SCET studies of low-p^ resummation and with the standard QCD approach. In the remain- 
ing sections we present fixed order results for the various perturbative functions appearing 
in the factorization theorem, anomalous dimension computations and RG equations, various 
internal consistency checks of the factorization theorem, and checks of the factorization by 
comparing with known fixed order QCD results. 

To recap, the factorization theorem for the Higgs transverse momentum p^ and rapidity 
Y has the form 

d 2 a 7r 2 f 1 dxi f 1 dx 2 f 1 dx[ f 1 dx' 2 

dp 2 T dY = 4(iV2 - i)2Q2 J ^77 J xi ^r J x2 ^r 

x E{x x x 2 Q 2 , fiQ- fi T )G lJ (a?i, x[, x 2 , x 2 , p T , Y, /i T )/ i/P (x / 1 , Hr)fj/p(x' 2 , fPp), 

(54) 
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where H and Q % i are perturbative coefficients and f\/p are the standard QCD proton PDFs. 
The indices i,j run over the gluon, quark, and anti-quark species. In our initial analysis we 
have ignored contributions from quark and anti-quark initiated Higgs production since they 
are numerically small, but it is straightforward to generalize the analysis and include these 
effects if desired. 

The factorization theorem achieves a clear separation of the physics associated with 
the scales Q ~ m h ^> p T ^> Aq C d which are now encoded in the hard function 
H(xiX2Q 2 , fiq; Ht), the perturbative function Q l ^(xi, x[, £2, x' 2 ,Pt, Y, fir), and the non- 
perturbative PDFs respectively. The hard coefficient H(xiX2Q 2 , Hq\ fir) is obtained by 
matching rif = 5 QCD onto SCET PT at the scale fiQ ~ Q ~ m h , as in Eqs. (27) and 
(42), followed by RG running down to the scale \lt ~ Pt- The RG running sums up the 
logarithms of rrih/pT- The perturbative function Q^{x\, x[,x 2 , x' 2 , px, Y, n T ) is a convolution 
of two perturbative collinear functions and an iSF(inverse Soft Function). The result for 
involves logarithms of Pt/^t and it is thus natural to choose fix ~ Pt to minimize these 
logarithms. The QCD PDFs fi/p(x,fiT) are also evaluated at this scale \xt ~ Pt and the 
logarithms of Aqcd/pt are summed via the standard DGLAP evolution of the PDFs. This 
picture is roughly summarized in Fig. [T] We note that in principle, the iSF could depend 
on an additional scale [is- In this analysis we set us = (J>t, since we see no evidence of large 
logarithms associated with the additional scale Q that must be resummed via an additional 
renormalization-group evolution in [is- We show this when calculating the iSF in Section [V] 
If a higher-order analysis indicates such a need for a separate fis, it is straightforward to 
include. 



The perturbative function Q lJ defined in Eq. (51) has an intricate structure on which we 
now elaborate. As seen in Eq. (51), is a convolution, in light-cone momenta and an 
impact parameter, over a n-collinear function X^7 g i7 n-collinear function and the iSF 
iS _1 . The l-n^n-gi functions are obtained from matching the iBFs (impact-parameter Beam 
Functions) B^f n onto the standard QCD PDFs as defined in Eq. (43) and are related to the 
iBFs as in Eq. (46). The iBFs are defined without a soft zero-bin subtraction but with an 
ultrasoft zero-bin subtraction. This is a consequence of the equivalence of soft zero-bin and 
soft subtractions [73TT75] which we have discussed in some detail in Section III C If one works 
with iBFs defined with a soft zero-bin subtraction, the purely collinear iBFs defined earlier 
in Section HI C then these objects would naturally match onto purely collinear PDFs which 
differ from the standard QCD PDFs by soft zero-bin subtraction terms. For more details 
on the relationship between standard PDFs and purely collinear PDFs see [7U [75] . 

The n-collinear and n-collinear iBFs describe the physics of collinear emissions in the n 
and n directions respectively but in addition, as already mentioned, they each encode the 
same soft region due to the absence of a soft zero-bin subtraction. As a result the soft 
region, which encodes the physics of soft emissions and soft virtual effects, ends up being 
double-counted between the two iBFs. The iSF plays the role of precisely subtracting a 
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soft region so that in the convolution over the two iBFs and the iSF there is only one soft 
region as required. After the iBFs are matched onto PDFs as in Eq. (43), we end up with 
a convolution over the two collinear functions 1^f n . gi anc ^ the iSF which defines as in 
Eq. ([5l|. 

The iBFs B n n are defined in covariant or non-singular gauges where the gauge poten- 
tial vanishes at infinity. Under non-singular gauge transformations, the iBFs are invariant. 
However, the iBFs are not invariant, as defined, under singular gauge transformations with 
boundary conditions where the gauge potential is non-vanishing at infinity. In singular 
gauges like the light-cone gauge, the definition of the iBF must be modified [4*4"144"T] with a 
transverse gauge link in order to be invariant under singular gauge transformations. How- 
ever, the restriction of the iBFs to covariant gauges suffices since the quantity which is 
a convolution over the two iBFs and the iSF is fully gauge invariant. This can be seen in 
two ways. First, since the iBFs and the iSF were obtained after applying the soft-collinear 
decoupling property of the SCET PT Lagrangian on the manifestly gauge invariant operator 
O defined in Eq. (21 ), the total convolution over the iBFs and the iSF must again be gauge 
invariant. Second, the product of the hard coefficient H and the perturbative function Q l \ 
appearing in the factorization theorem in Eq. (54), just gives the perturbative cross-section 
for Higgs + multi-parton production as is manifest before applying the soft-collinear decou- 
pling transformation. Together with the gauge invariance of H, the quantity must again 
be gauge invariant. The additional functional dependence of the iBF on the variable t n also 
removes the singularity that occurs when this variable is fully integrated over, as it is in the 
transverse-momentum dependent PDF. In our factorization theorem this variable is instead 
constrained by the external kinematics of the process. 

The factorization of Q 1 ^ into the convolution X^ g i ® X n . g j a/3 <S> can be quite useful 
in the implementation of higher order radiative corrections in oi s (pt)- This is because the 
computation of the individual functions 1^. gi , defined in covariant gauges, and S~ l is much 
simpler than the full computation of Q % i without factorization. 

Another property of the iBFs B^n an d the matching coefficients X^ n . gi is their index 
structure. The indices a, /3 run over the transverse directions. Note that these indices 
remain uncontracted in each iBF due to their non-trivial dependence of the iBFs on the 
impact parameter 6". Instead, these indices are contracted between the n-collinear and the 
n-collinear functions as seen in Eqs. (41) and (51). This contraction of the indices between 
the n and n collinear functions is necessary to reproduce dot products between the transverse 
momenta of final state particles in the n and n directions which will appear at higher orders 
in perturbation theory. 

We next discuss the formulation of factorization in impact-parameter space vs momentum 
space. So, far we have discussed the factorization theorem in terms of the impact-parameter 
space iBFs and iSF as seen in Eq. (51). In particular, the matching of the iBFs onto the 
standard QCD PDFs, as defined in Eq. (43), must be done in impact-parameter space. It is 
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only in impact-parameter space that the infrared singularities of the iBFs become manifest 
as poles in dimensional regularization, and can therefore be seen to be the same as those 
of the standard PDFs. However, once the impact-parameter space matching coefficients 



are obtained one can Fourier transform to momentum space as in Eq. (52) and obtain the 



perturbative function entirely in terms of the momentum space matching coefficients and 



iSF as seen in Eq. (53). In this momentum space factorization there is no reference to any 
impact-parameter. This form also makes manifest that the Landau poles which appear in 
the impact-parameter space CSS formulation [31] are avoided here. We discuss this point 
further below in our comparison with the standard approach. 



B. Comparison with CSS approach 

The classic QCD analysis of resummation of low transverse momentum logarithms ex- 
presses the cross section in the low px region as 



dpx dY 

x exp 




d\i 



a.b 

2 



\n%A{a s {^)) + B{a s {^)) 



(55) 



The sum is over parton species labeled by a, b, while xa,b denote the equivalent parton 
fractions xa,b = e ±Y rrih/Q respectively. The functions A, B, and C have perturbative 
expansions in a s , while b is an arbitrary constant chosen for computational convenience. 

One significant difference between this result and our approach outlined in the previous 
section is the appearance of the Landau pole of the strong coupling constant when fi 2 = 
in the exponent. To deal with this singularity, several modifications of this formula are 
employed, including a deformation of the b± integration contour [361 EE [77] and the intro- 
duction of a phenomenological model to cut off the b± — > oo region [78]. In our approach 
the most natural choice for the scale which controls the lower limit of the RG evolution 
is hl = pt- This can also be understood by noting that the perturbative function is 
independent of the impact parameter, in both the impact-parameter and momentum-space 
formulations of the factorization theorem, and depends on pt and \x and no other dimen- 
sioned scales. Furthermore from the structure of the factorization theorem, we see that the 
logarithms of ruh/pT are summed by the RG evolution of the hard coefficient H(xi, X2Q 2 , /i) 
which multiplies the function Q l i and also has no reference to an impact parameter. In 
the effective theory, non-perturbative effects such as those indicated by the appearance of 
the Landau pole are encoded in operators suppressed by Aq C d/pt- When p T ~ Aq C d, 
the expansion in this parameter breaks down, and a model of Q 1 ^ fit to data can be used 
analogous to the standard approach. However, no reference to a non-pertubative function 
is needed above Aq C d- Previous comparisons of 6-space and momentums-space resumma- 
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tion formalisms have indicated numerical agreement between the obtained results down to 
Pt ~ few GeV [79J. At this stage, power-suppressed operators presumably give important 
contributions. The use of SCET allows such effects to be studied in a systematic way. The 
avoidance of the Landau singularity also simplifies the matching of the resummed result to 
the fixed-order expression. In the usual approach, a large cancellation between the resummed 
component and the fixed-order QCD contribution occurs, leading to potential instabilities 
in the matched distribution. This cancellation typically occurs numerically because of the 
introduction of a non-perturbative model for the large b± region. Since it can be arranged 
analytically if the b± integrals can be done, avoidance of the Landau pole is useful for this 
purpose also (we note that the matching to fixed-order QCD results can be made smoother 
in the CSS approach by an appropriate change in the logarithm resummed [35]). We note 
that because of the effective theory expansion in Aq C d/pt, it is not possible to calculate 
exactly the pt = result for very large Q as in the 6-space approach [28J. However, the 
model of the non-perturbative region can be arranged to reproduce this result [79J . 



C. Comparison with previous SCET analyses 



Resummation of the of the Higgs production transverse momentum distribution has been 
studied previously in [13] by a perturbative matching calculation onto SCET followed RG 
running. An explicit all orders factorization formula at the pt scale, corresponding to 
our Eq.Q, was written in [12]. The structure of our factorization theorem differs from 
that in [42J in several key ways. Our factorization theorem involves a convolution over 
the residual light-cone momenta that connects the collinear and soft sectors necessary by 
momentum conservation which is absent in [42J. This additional dependence is required 
in order to extend the factorization and resummation to more exclusive quantities where 
final-state hadronic quantities are also restricted, as is clear from the dependence on the 
light-cone momenta in the beam functions of Ref. [41J. The iBF in our factorization theorem 
has transverse indices which are contracted between the n-collinear and ra-collinear sectors 
necessary to generate dot products between the transverse momenta of partons in the n and 
n directions in the final state. This index structure is not present in [12]. We disagree with 
this result, as it is not the structure that arises from the SCET analysis. In our formulation, 
the matching onto the standard QCD PDFs was done separately in the n and n sectors 
making full use of the SCET factorization. The work of [12] does the matching to PDFs by 
matching the full product of soft and collinear functions onto the QCD result, as required 
if the index structure of the collinear sectors is dropped. Finally, our work incorporates the 
formalism of zero-bin subtractions necessary to avoid double counting the soft region which 
was not addressed in Ref. [42J. 
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FIG. 3: The SCET diagrams contributing to the calculation of the Wilson coefficient C{to\ , uj 2 , /x) . 
The purple cross denotes the n and n collinear Wilson lines and the soft Wilson lines, while gluons 
with lines drawn through them are collinear gluons as the n and n labels indicate. The S label 
denotes a soft gluon in the first diagram. 



V. FIXED ORDER ANALYSIS 



In this section we give details of the perturbative computation of the hard Wilson coef- 
ficient, the n-collinear and n-collinear iBFs, and the iSF which appear in the factorization 
theorem. We emphasize that the iBF and iSF derived here are universal, and can be applied 
to other gluon-initiated processes at hadron colliders. 



A. Calculation of the QCD -> SCET PT Wilson coefficient 

We begin by discussing the matching of QCD onto SCET PT in order to extract the Wilson 
coefficient C(ui, u 2 , fj)- The Wilson coefficient can be extracted from the relation presented 



in Eq. (27) by computing radiative corrections to the matrix elements of both the QCD and 
SCET operators and encode their difference in C(ui,u 2 ). For the tree level and one loop 
matching one can compute the matrix elements {h\OQCD\Pi,P2) and (h\0\pi,p 2 ) in QCD 
and SCET PT respectively where p\^ = n ■ P\\ and p 2 = n ■ Vt^ denote the momenta of the 
initial state n-collinear and ra-collinear gluons. The diagrams contributing at next-to-leading 
order in a s in SCET PT are shown in Fig. [3j Labeling these graphs from left to right they 
take the form 

Fig. [3}i = V a (p 1 ,p 2 ) 0{n -pi,n -p 2 ), 

Fig. [3)3 = [V b {p 1 )-V w ip 1 )}0(n-p 1 ,n-p 2 ), 

Fig. [3^ = [V b (p 2 )-V b0 (p2)]O(n-p 1 ,n-p 2 ), (56) 
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+ c.c 



FIG. 4: The diagrams contributing to the next-to-leading order jet function. The purple cross 
denotes the collinear Wilson lines associated with the B± field. The blob in the left-most diagram 
denotes the wave- function graphs. We note that the momentum p\ is incoming on the left-hand 
side of the cut and outgoing on the right. 



so that the SCET PT operator is multiplicatively renormalized. With on-shell external gluons 
and using Feynman gauge, the quantities V a ,b,bo take the form 



V b (pi) 
Vmipx) 



-ig 2 c A ) 



-ig 2 c A ) 



:-*g 2 c A ) 



d d 



(27r) d £ 2 n-£ n-£ 
d d £ (n 



2 + (n ■ pi) 2 + n ■ 



n ■ pi 



(27r) d 
d d i 



£ 2 (£ + Px) 2 n- (i + Pi)n-£ 
2 



(2n) d £ 2 n-£ n ■ £' 



(57) 



We note that the collinear graphs in Figs. [3p and |3p require a zero-bin [72] subtraction 
given by the V^o term in order to avoid over-counting the soft region. These integrals are all 
scaleless and vanish in dimensional regularization. The Wilson coefficient can therefore be 
extracted directly from the finite part of the full QCD result in dimensional regularization 
and is well known in the literature [151 80J . Through next-to-leading order, it is given by 



C(n -p x n ■ p 2 ,fi) 



cn ■ pin ■ P2 

v 



Ti- 



ll + — - ln z 

6 



n ■ pin ■ P2 

V 2 



(58) 



where the first term is just the result of tree level matching quoted earlier in Eq. (28) and 



we have used the property C(co x , U2, /x) = C(uiU2, /i) to write the LHS above. We note that 
the 11/2 in the next-to-leading order expression arises from integrating out the top quark 
loop to produce an effective Higgs-glue vertex. The hard Wilson coefficient H(xiX 2 Q 2 , /i) = 
\C(xiX2Q 2 , /i)| 2 appearing in the factorization theorem can be obtained from the above 
equation at next-to-leading order. 
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B. Calculation of the iBF 



In this section we present results for the calculation of the iBFs as defined in Eqs. (36), 



(33), and (42). We compute the n-collinear iBF by inserting a complete set of states as 

initial pols. X n 

x (X n \5(P-x 1 n- P i)gBt ±a (0)]\px), 



(59) 

and then computing the product of matrix elements. Recall that denotes the iBF 
without a soft zero-bin subtraction as opposed to which is defined with a soft zero-bin 
subtraction. An analogous expression exists for the n-collinear iBF. In this section we focus 
on the n-collinear iBF, since the n-collinear iBF can be calculated in an analogous fashion. 
The lowest order result for the iBF is obtained by choosing \X n ) = |0) and computing the 
tree level matrix elements to get 

B^ixxXM = B^(x u t+,b ± ) = -g 2 gf5(t+)5(l - Xl ). (60) 

Higher order radiative corrections to the term in the iBF with \X n ) = |0) correspond to pure 
virtual corrections with no real emissions in the final state. The one loop virtual correction 
corresponds to the first diagram in Fig. [4] and its conjugate and is given by 

B^( Xl ,t+,b ± ) = B^( Xl ,4,b ± )\v b ( Pl ) + V b (- Pl )} +w.f., (61) 



where V b is given in Eq. (57) and w.f. denotes the external wave- function graphs. Thus, 
this virtual correction is also scaleless and vanishes in pure dimensional regularization. 

Next we compute the single real gluon emission contribution to the iBF corresponding to 
a single gluon with momentum p g in the final state: \X n ) = \p g ). The single gluon emission 
contribution to the iBF corresponds to the last three diagrams in Fig. |1J Only the second 
diagram in Fig. [4] contributes if we use the physical gluon polarization sum 

Y,^V(P 9 ) = ST + 9 - 9 ■ (62) 

pols. " 

Direct calculation of this graph gives 

B^(x 1: t+,b ± ) = [ d d p g 8(pl)5(pT-x 1 Q-p;)5( t i-p+)e^ 




(p ^ 2 



J 9± 



X l Pg± J 



(63) 
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where we work with the convention d = 4 — 2e. After performing the integrations over the 
light cone momenta p^ and p~ the single gluon emission contribution to the iBF can be 
brought into the form 

2g A ^C A 



x \gf 



27r) d - 1 
1 -xi -+ 



1 



1 — X\ 



(d-2 



1 ~ Xl Pg±Pg± 

x\ p 2 g± 



(64) 



Three-dimensional quantities have been denoted by a vector symbol. The result is the 
same in any gauge due to the structure of the Wilson-line coupling. We note that if 
was integrated over, as it is when considering transverse-momentum dependent PDFs, 
singularities would occur at the integration boundaries. This is avoided in the iBF. The 
variable t£ is instead set by the external kinematics through the convolution structure of 
the factorization theorem, as we will see in a later section when comparing to fixed-order 
QCD. 

In order to perform the integral appearing in Eq. (64), we require the following two 
integrals: 



S[tt(l-x 1 )-p'] . 



t+ 



5[t+(l - Xl ) - p 2 ± ] ^ ± , fg± (b ± ■ p g± ) 



t+ 



-(2 7 r) 1 - e &l(t+)- 1 - £ / 2 (l - Xl r' 2 J_ c [b ±y /t+(l-x 1 )] - (d - 3) 



t+(l-x 1 ) 



(65) 



and we have used the notation = b±. To utilize these results, it is most useful to expand 
the integral over the second piece of Eq. (64) in terms of form factors. We can schematically 
write this as 



d 2 ~ 2 %/-^ = g 1 b\ g f + g 2 bib{. 



(66) 



We can express these in terms of the two integrals calculated in Eq. (65): 

_ gZi - X 2 _ bjl, -{d- 2)X 2 
^-&i(d-3)' 92 ~ 6*(d-3) ■ (67) 

Furthermore it is convenient to parameterize the iBF in terms of two form factors T\ and 
J-2 as 



B^(x x ,t+, b ± , t i) = T 1 (x 1 ,t+, b^^gf + F 2 (x 1 ,t+, b ±lt j.) 



gf + (d-2)- 



(68) 
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Note that if the indices a, (3 are contracted the term proportional to Ti vanishes. Using 
Eqs. (64), (65), (66), and (67) we find the following results for the form factors T\ and Ti 



2g^ 2e C A 

_ (27T) 4 *- 1 

1 - Xi H t, h 



l(2ny-^ ± (t+)- 1 -^(l - Xl )-^J- e [b ± y/t+(l- Xl )] 



x\ 



1 — X\ 



2g A ^C A 
(2n) d - 1 



-(2 7 r) 1 - e 61(t+)- 1 - e / 2 (l - Xl )-^J 2 . e [b ± y/t+(l- Xl )] 



1 — X\ 



(69) 



The form factor T\ has ultraviolet and infrared divergences that are regulated by e in our 
pure dimensional regularization calculation. On the other hand, the form factor T% is finite. 
It is convenient to remove the fractional powers of the Bessel-function arguments by writing 
them in terms of hypergeometric functions in the following way: 



Uz) 



0*1 1 + v\ 



4 



(70) 



Since qFi{v, 0) = 1, and therefore J%{z) ~ z 2 , it is clear that T<i is finite and that we can 
safely set e = to get 



57T X\ 



bitm-xt) 



(71) 



For the form factor J-'i, we first rewrite its expression in Eq. (|69j) as 



Q 5 



i /t+ 



x o*i 1 - e; 



7T e 

,2 



Q 2 \Q 2 



-l-e 



1 - X! 1 

1 — Xi H , f 



1 — Xi 



ftjC(i-si) 

4 



(72) 



where we have have switched to an MS" definition of \x by replacing /i 2e — >■ /i 2e (47r)~ e e <E7 . 
We have also used the partonic center-of-mass energy squared Q 2 = X\X 2 Q 2 to form dimen- 
sionless ratios where required; this is convenient for the comparison with fixed-order QCD 
calculations, which is discussed in a later section. We note that the ratio of scales which 
appears in this expression is t^/fi 2 , indicating that we should choose \x ~ pr to minimize 
logarithms since in the partonic calculation ~ p\. We write the expansion of the form 
factor T\ in e as 



+ ^1;C 



(73) 
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and derive the following expressions 



Fx* 



.^C7^(t+) < 5(l-x 1 ), 

7T 



a. 



= ^C^- <^ 5(1 - X a 
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Q 2 



5(t; 
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■1 1 
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-5(4)5(1 - X!^ 
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-5(t+) 5(1 -x x ) + 



ti- 



ln 



Xi + 



.7.'! 



XT 
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4 
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+ -5(t n )5(l -xi)ln 
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x\ 
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1 — X 
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oFl ( 1; _SSt(iz-)) 



Xi + 



X\ 



XT 



1 - X! 



In 



+ 5(1 - xi) 



(74) 



The result for the single gluon emission contribution to the n-collinear iBF in pure dimen- 
sional regularization is given by Eqs. (68), (71), (73), and (74). The expression for the 
n-collinear iBF is obtained by repacing 4 — >* an d x i ~~ >* x 2 i n the expression for the 
n-collinear iBF. 



C. Soft function: real emission 



In this section we give results for the computation of the iSF which was defined earlier as 



S (ui,u 2 ,b±,fi) 



db+db ih+Gjl 
16tt 2 



/V 6 ** /2 s- l {b + XM, 



(75) 



where the position space soft function that appears on the RHS above is defined in Eq. (33) 



The iSF in the factorization theorem, before doing the Higgs phase space integrals, has the 
arguments 

Wx-Ph- k ni u 2 = uj 2 -Ph- k+, (76) 



as seen in Eq. (47). For convenience we introduce the notation 

= Q(U2- P +), C 



t: 



Qh$, t- = Qk~, tr x = Q(u; 2 - p+), tr x = Q(ux - Ph), (77) 

which we will often use in this section. We compute the iSF by inserting a complete set of 
soft states in the position space soft function as 

5(6, /x) = ^(0|f [Tr (s n T D YiS n T c Sl) (b)} \X S )(X S \T [ir (s n T c SlS n T D Sl) (0)1 |0), 



(78) 
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FIG. 5: Example diagrams contributing to the next-to- leading order iSF. The four lines at each 
vertex schematically denote the soft Wilson lines appearing in the definition of the iSF The 
diagram on the left corresponds to a virtual correction to the iSF and the diagram on the right 
corresponds to a real emission as seen by the cut through the gluon. 



and compute the product of matrix elements and use these results in Eq. (75). Through 



next-to-leading order in the QCD coupling, the position space inverse soft function S 1 (b) 



that appears in Eq. (75) is obtained by inserting an overall minus sign in the 0(a s ) correction 



to the soft function S(b) of Eq. (78). 



The lowest order result for the iSF comes from choosing \X S ) = |0) and computing the 
tree level result which gives 

N 2 - 1 

5 -Ko) (C - _ C * _ £ - 6±> M) = ^—Q 2 6 (tr x - t- n )5{tT x - t+). (79) 

Higher order corrections to the term with \X S ) = |0) corresponds to virtual graphs with no 
real emissions in the final state. At one loop, the virtual corrections correspond to the first 
diagram and its permutations in Fig. [5] which gives the result 

5- lv «(wi-Ph,^-^,6j.) = -5 v «(w 1 -pfc,w ! ,-p+ b ± ) 

= - p -,u 2 - p+,b ± )(-2tg 2 C A )I s 

(80) 



where I s is the scaleless integral 

d d l 



(81) 



(2vr) d (£ 2 +i0)(n-£- iO) (n-£ + i0)' 

and vanishes in pure dimensional regularization. 

Next we compute the contribution to the iSF from the real emission of an soft gluon 
corresponding to choosing \X S ) = \k) for a gluon of momentum k, as shown in the second 
diagram of Fig. |5j Explicit computation gives 

N 2 c -lg 2 ^C A f db + db- eib+ a l/2eib -z 2/2 f d d k5 ^ 4 e _ ib . k 



(27r) d - 1 J 16tt 2 J v J k+k~ 

(82) 
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Switching to an MS definition of fi and performing integrals as before, we can derive the 
following expression: 



N 2 - 1 a s C A 



■„/■■ (l r. 1,1 



-a uj 1 u 2 



4 7T r(i-e)^ 1 z u x v 4 y 

The expansion in e proceeds identically to that for the iBF. Defining the expansion 
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we arrive at the result 
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We have used the scale Q to define dimensionless ratios in logarithms and plus distributions, 
as in the calculation of the iBF. We note that when this result is plugged into the full 
expression for Q l i in Eq. (49), the logarithm that appears is ln iu/(n 2 Q 2 ) , where t,u are 
the partonic Mandelstam variables. This indicates that we should choose \i ~ pr in the soft 



function, and is the motivation for our decision to set \is = in Section [TV 

D. iBFs to PDFs 

In order to be able to sensibly match the iBFs onto the PDFs with a finite Wilson 



coefficient, as indicated in Eq. (43), the infrared divergences in the iBF and PDF must 
match. In this section, we show that this can be manifestly seen from the integrand level 
expressions for the iBF and PDF at the order we are working. We first note that the tree 
level expression for the iBF to PDF matching is given by 



(86) 
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which can be straightforwardly verified from Eq. (43). Next we consider the integrand 



level expression for the real gluon emission contribution to the iBF in Eq. (64) where the 



integrals over the light cone momentum components of the final state gluon have already 
been performed. The infrared divergences now occur in the limit where the perpendicular 
momentum of the final state gluon vanishes p g ± — > 0. In this limit, the iBF real-emission 



contribution in Eq. (|64|) becomes 



-gfs(t; 



IR 



(2tt; 



d-1 



1 — X\ 



1 — X\ 



x\ 



(87) 

A calculation of the real-emission contribution to the PDF defined in Eq. (44) can be done 



inserting a single-gluon state between the field operators in the same way as was done for 
the iBF. The result at the integrand level is 



fifpM 



2g 2 ^C A fd 2 - 2 % ± 
(2< 



d-l 



X\ + 



1 



+ 



from which we see that 



[fW( Il) C l ii,4 = -S , « 



Xl 



(88) 



(89) 



Thus, the infrared divergences which appear upon integration over p g ± are identical for 
both the iBF and the PDF, up to an overall multiplicative factor. Similarly, we note from 



Eqs. (61) and (60) that the contribution of one loop virtual corrections to the iBF take the 
form 



Xl 



(90) 



so that once again the infrared structure is the same for the iBF and PDF virtual contri- 
butions up to the same multiplicative factor. Thus, the infrared divergences in the sum of 



the real and virtual graphs for the iBF can be matched to those in the PDF as in Eq. (43) 



with the finite Wilson coefficient given in Eq. (86). Additional contributions to the Wilson 



coefficient will come from the finite terms in the iBF. 

In fact, since the PDFs are scaleless and we have shown that the infrared divergences 
of the iBF and the PDF match up, the Wilson coefficient for the iBF to PDF matching is 
just given by the finite part of the iBF in pure dimensional regularization as indicated in 



Eq. (46). For a more detailed explanation of this point see [81J. Using Eqs. (46) and (68) 



we can now write the iBF to PDF matching coefficient at the first order beyond tree level 
as 



x x 



AO = ~ x i { ^i;o(^r s x'it+ , b ± , fj)g 



X, 



+ T 2 (^-,x^ 



i i ^l^n ' 



Xi 



gf + (d-2)- 



(91) 
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where the form factors J-i-o and T 2 are given in Eq. (74) and (71) respectively. We perform 



an explicit check of this expression in section fVTT] by using this expression in the factorization 
theorem and comparing the fixed order result in full QCD. 



VI. RUNNING 



In the factorization formula of Eq. (54 ), the logarithms of rrih/pT are summed via the RG 



evolution of the hard coefficient H(xix 2 Q 2 , //) between the scales /jq ~ and /jt ~ Pt- 
The factor H(xix 2 Q 2 , Hq] (J>t) m the factorization formula denotes this RG evolved hard 
coefficient. The logarithms of ^.qcd/pt are summed via the DGLAP evolution of the PDFs 
up to the \it ~ Pt scale. The DGLAP evolved PDFs are denoted by fi/p(x, //t)- In this 
section we give details of this RG running above and below the px scale. 

A. Running above pt 



The anomalous dimension of the hard coefficient H can be obtained from the anomalous 
dimension of the operator 



which was first defined in Eq. (21). The RG equation for O has the form 



d = _ 7o o, 

an 



where the anomalous dimension 70 is given by 

7o 



1 d 
—fx—Zo, 
Zjq ail 



(92) 



(93) 



(94) 



and Zq is the renormalization constant that relates the bare operator O b to the renormalized 
operator O by 



O b = Z 0. 



(95) 



We use the notation u\ = x\n-p\ = n-pi, 002 = X2n-p2 = n-p 2 and n-p\n-pi = X1X2Q 2 = Q 2 . 
The anomalous dimension of 0(n ■ p\,n ■ p 2 ) is typically written in the general form 



-Q2 

70 = r [a a ] In — — + 7o[a s ], 
where the cusp (To) and non-cusp (70) pieces have expansions in a s as 
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Equivalently, one can think of the running of the Wilson coefficient C of the operator O 
which has the RG equation 



l*-j~ C = 7 C C, 



1 d 7 <7 - 7-1 



7c = lo- 



rn 



The running of absolute value squared of the Wilson coefficient H(Q 2 ,fi) = \C(Q 2 , /i)| 2 , 
which is the quantity that appears in the factorization theorem, is given by 

^H = lH H ) lH = -A_ fi A-Z H , Z H = (Z + Z* )- 1 7H = 7 c + 7c *, (99) 



so that 



Q 2 

1h = T H [a 8 ] In — + IhI^s 



where V h and 7^ have the a s expansions 
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IhWs] = ^7o J + 
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47rJ 



(100) 



(101) 



The result for the anomalous dimension 70 is well-known [TBI E31 IHDl E21 IH3]- Since the 
SCET PT graphs that renormalize the operator O are scaleless, as discussed in Section |Vj the 
anomalous dimension can be read off [81] from the infrared pole structure of loop correc- 



tions [80j [82] to the QCD operator Oqcd defined in Eq. (20). In this section we reproduce 
this known result at one loop by explicit calculation of the SCET PT graphs using off-shell 
external gluons to regulate the infrared divergences and dimensional regularization to reg- 
ulate the ultraviolet divergences. We do this calculation in both Feynman gauge and using 
the background field method as an additional check on our calculation. 



1. Anomalous dimension computation in Feynman gauge 



At one loop, there are three graphs in SCET PT , in addition to the wave function graphs, 
which determine the anomalous dimension 70 = 7 C , as shown in Fig. [3j The first graph 
involves a soft gluon loop emanating from the soft Wilson lines and the remaining two 
graphs involve a (n,n)-collinear gluon originating from the (n,n)-collinear Wilson line and 
attaching to the (n,n)-collinear gluon. These diagrams are given by 

Fig. [3^ = V a (pi,p 2 )0(n-p 1 ,n-p 2 ), 

Fig. [3b = [V b (p 1 )-V b0 {p 1 )]O{n-p 1 ,n-p 2 ), 

Fig. [3b = [V b (p 2 )-V b0 {p 2 )]O{n-p 1 ,n-p 2 ), 



(102) 
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where the quantities V a (pi,p2), Vb(pi), and V&o(pi) with off-shell momenta pi^ are given by 

d d £ 2n ■ pin ■ p 2 



V a (pi,p 2 ) 

Ho(pi) 



(ig 2 C A ) 



-ig 2 c A ) 



-ig 2 c A ) 



(2n) d (£ 2 + z'O) [n ■ £ n ■ p l + p\ + £0] [-n ■ £ n ■ p 2 + pi + iO] 
d d £ (n ■ £) 2 + (n • pi) 2 + n ■ £n ■ pi 

iO] [n-£ + iO] 



(2n) d (£ 2 + z'O) [(£ + pi) 2 + zO] [n ■ (£ + 
eft 2n • pi 

(27r) d (£ 2 + i0) (n- £n-pi+p 2 + i0){n- £ + iO] 
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and note that in the limit p 2 ,pl — > these integrals reduce to the expected results of Eq.(57). 
Note that the (n,fz)-collinear graphs require the zero-bin [72] subtraction terms Ho(Pi) an d 
Hofe) respectively. These zero-bin subtractions are necessary to remove the soft region in 
the collinear graphs and thus avoid double counting the soft region. The results for these 
graphs and the zero-bin subtraction term are 
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n ■ p\J J J 
We note that the off- 



where we have used an off-shell infrared regulator A 2 = p\ 
shellness of the external gluons does not completely regulate the infrared divergences indi- 
vidually in Vb and V&o- Furthermore, V b contains a UV-IR mixed double pole. However, in 
the zero-bin subtracted combination V b — V&o the infrared divergences are completely regu- 
lated by off-shell external gluons and all pole terms are pure UV. This is critical for being 
able to properly extract the anomalous dimension from the pole structure. The sum of these 
graphs is 



Fig. Hi + Fig. [3b + Fig. 03 
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Note that the mixed UV-IR terms proportional to — In A 2 cancel between the soft and 
collinear graphs so that the UV and IR divergences are clearly separated in the sum of all 
graphs. Along with the gluon wave-function and strong coupling renormalization constants, 
we can now extract the anomalous dimension of O. The bare operator O b obtained by 
expanding the bare B\ fields in the strong coupling and keeping the leading term, can be 
related to the renormalized operator O as 



6 b = 9flu hTr [(gA^UigA^U] =Z Q 6 



(106) 
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so that the renormalized operator can be written as 

Z 2 Z 

6 = g^hTr [(gA^UigA^U] + - lj^Tr [(gA^^A^)^], 

= g^hTr [(gA» ± U(gA» ± U} + (25Z g + 5Z, - 5Z )g^hTr [(^UG^JJ +■■■ , 

(107) 

where Z3 and Z g are the gluon wave function renormalization and coupling renormalization 
constants respectively 

A b = Z 3 A, g b = gZ gf i e , 



and we have defined 
Z 3 



SZ, 



l + SZ n 



1 + 5Z 
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At one loop, in the MS scheme, the gluon wave function and the strong coupling renormal- 
ization constants are known 
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We see from the above results that in the MS scheme 5Zo at one loop is given by 

8Z$ = 2 SZf + 6Z® + [K(pi,p 2 ) + V h $ x ) - V m {p x ) + V b (p 2 ) - V b0 (p 2 )] d . v , (111) 

where the subscript "div" above indicates the UV divergent part of the sum of graphs. 
Putting these results together we get the one loop result 
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Using N c = Ca = 3 this result simplifies to 
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from which we get the well known result for the one loop anomalous dimension 7^ 
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In the notation of Eqs. (96) and (97) this result for the one loop anomalous dimension 
corresponds to 



F° - 4 



7o° = 0. 



(115) 



Thus, the non-cusp piece 7^ vanishes. Putting this together with the notation in 
Eqs. (99), (100), and (100) the anomalous dimension of the hard coefficient H at one loop is 
by 



= -2r ° = -8, 
and also contains only a cusp contribution. 



7o" 



-27o° = 0, 



(116) 
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2. Anomalous dimension computation in background field method 

One can also compute the anomalous dimension using the background field method with 
Feynman gauge for internal lines in the loop graphs, which we show as an additional check 
on our SCET computational techniques. In this method, the combination gA^ ± is not 
renormalized. In this case, the result for the soft graph V a g is unaffected 
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and the zero-bin subtraction term is unaffected 

n b o fg = Ho. (H9) 

One can now extract the anomalous dimension following the same procedure as in the 
Feynman gauge computation of the last section. The bare operator obtained at leading 
order by expanding the bare B b ± field can again be written in terms of the renormalized 



operator as in Eq. (106), so that the renormalized operator can be written as 

1 



O 



g^hTr [{gA^J^gA^)^] + - ljg^hTr [(gA% ± ) wl (gA^ 
g^hTr [Ki)«Msi) u J - 6Z g^hTr [(gA^^gA^)^] + 



(120) 



Note that in the background field method, compared to Eq. (107), the wave function (Z3) 



and coupling constant (Z g ) renormalization constants cancel out and do not appear. In 
other words, the combination gA^ n± does not undergo renormalization in the background 
field gauge. Thus, we get 



SZ. 



(i) 
o 



Puh) + v b hts (Pi) - HTOO + H btg (p 2 ) - Ho s (p 2 )] 



div' 
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which gives the same result as in Eq. (113) and thus leads to the same anomalous dimension 



as derived from the Feynman gauge calculation. 



B. Running below pt 



Logarithms of ^.qcd/pt are summed via the standard DGLAP equations for the PDFs. 
For completeness we write down the RG equation for the gluon PDF at leading order 
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where P g ^ g (z) is the standard gluon splitting function and we have ignored contributions 
from the quark and antiquark PDFs. This DGLAP equation is used to evolve the gluon 
PDF up to the fix ~ Pt scale. These evolved PDFs, denoted by f g /p(x, /xt), are the ones 



that appear in the factorization theorem as in Eq. (54). 



VII. CONSISTENCY CHECKS 
A. Cross-section comparisons with full QCD 

In this section we show that the iBFs and iSF computed in fixed-order perturbation 
theory in Section [V] reproduce the QCD cross section up to terms suppressed by the ratio 
PT/m h . 

1. Leading-order cross section 
We begin by reproducing the leading-order cross section of perturbative QCD for gg — > h, 



with no final state gluon emissions, using our factorization formula in Eq. (50). We utilize 



the leading-order expressions for the Wilson coefficient H, the iBFs and the iSF: 

'cxix 2 Q 2 ^ 2 



2^ 



H {0 \ Xl x 2 Q 

ln,gg ( — ,t n ,b±) = -g g/5 {t n )5 1 - — , 

x 2 dj 2 
-imax y.+ j-max -i— T\j2 i 

S~ m ( n n , n Q \ b ± ) = -^—Q 2 S(tT X -tn)S(tr X -tt). (123) 



Inserting these results into Eq. (50), the integrals over Xip, tn, tn, and b± can be easily 
performed to give 

d 2 a^ n /« s \2 r 

dudt* = 576v 2 \) / rfx i rfx 2/ 9 /p(a;i)/(,/p(^2)^(l-^)5(M-m^+Xi(5 2 )5(t-m^+X2(5 2 ). 

(124) 

We have introduced the variable z = m\j \x1X2Q 2 ) = ml/Q 2 , which measures the amount 
of energy released into final-state particles besides the Higgs (for z — 1, all energy in the 
partonic scattering process goes to the Higgs). It is convenient for this expression and for 
later results to switch to partonic Mandelstam variables u, t defined by u — m\ = X2{u — m 2 l ) 
and t — ml = X\(t — m\), and introduce a partonic differential cross section via 

dxi dx 2 f g /p(xi)f g /p(x 2 ) / dudi ——- ( 125 ) 

J du dt 
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The partonic cross section at leading order takes the form 

d 2 a^ 7T /a £ n 2 

7T 



(126) 



c£u dt 576f 2 

Upon integration over the variables u and t, this reproduces the leading order result of 



Eq. (13). 



2. Next-to-leading order cross section 

The next-to-leading order prediction of QCD has contributions from virtual corrections 
to gg — > h and from the emission of one final state gluon gg — > gh. To compare with this 



next-to-leading order prediction of QCD, we expand out all factors appearing in Eq. (50) 



through next-to-leading order in perturbation theory. After expansion, we can write the 
cross section in the schematic form 

(127) 

The next-to-leading order expressions for the terms appearing in this cross section can be 



du dt 



found in Eqs. (58), (85), and (91). Before explicitly computing the various contributions, 



we note several simplifying features that appear when the iSF and iBF terms are combined. 

• When the combinations Xn^X^ an d X^X^ are computed, a contraction of the form 
9°L-£a8 occurs ) since one of the two Wilson coefficients contributes only at leading 



order. This contraction removes the form factor Ti in Eq. (91), indicating that this 



structure contributes only at higher orders in perturbation theory. 

When the sum is performed, 

several terms that appear in the individual expressions cancel. In particular, it is 
simple to check that the following combination that appears in the the cross section 
vanishes upon integration over t n and tn'- 



x 




+ <H*n 



After making these simplifications, we have four contributions to the cross section at next- 
to-leading order. We now discuss each contribution in detail. To simplify the notaton, we 
will set fi 2 = Q 2 = x l x 2 Q 2 , where the arguments of the PDFs. 



41 



1. H^X^Xfj^ S This contribution has the same structure as the leading-order 



expression; using the Wilson coefficient in Eq. (58), we find the result 

7T 



a s \ 3 



11 7tt 2 



5(1 -z)5(u)5(t). 
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dud* 576v 2 V - ) '|2 12 
After integration over the Mandelstam variables, the contribution of this term to the 
total partonic cross section is clearly 

2 
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7T 



576w 2 
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2. HWZWZ^S-W: After implementing the simplifications described above, the only 



terms remaining from the iSF in Eq. (85 ) that contribute to the cross section are given 

by 
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The first term gives a contribution proportional to the leading order cross section. 



When computing the contribution of the second term to Eq. (50 ), the following identity 
is needed: 

bly\ 1 



db±b± J (k h ±b±) Fi 1 
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kh 



-5(k h± - y/y). 



(131) 



'0 V "* / 

This follows immediately from writing the hypergeometric function in terms of a Bessel 



function using Eq. (70), and then applying the orthogonality relation between J 
functions of different arguments. After performing the bj_ integral using the identity 
above, and expressing in terms of Mandelstam invariants using their definitions 
in Section |Vl we find the following contribution to the cross section: 
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To find the total cross section resulting from this term, we note that the range of the 
t integration is — Q 2 (l — z) < t < 0. The definition of the plus distribution l/[x]+ as 
the O(e ) term in the expansion of £~ 1-e allows us to derive the integral 
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Using this result, it is straightforward to derive 
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3. H^X^X^S- 1 ^: 



The terms from the next-to-leading order Wilson coefficient X n in Eq. (|91|) that con 

Q 2 



tribute to the cross section in Eq. ( |50| are 

1 



T (l,c)aP 
n,gg 



X F! 



x ln(l — xx/x-y) + 



-g C A — 9\ . 

7T I Q2 

J.2 (j-max 4.+\(4.max j.— ' 

u ±\ b n i n)\ i n 

4Q 2 

ln(l — Xi/x'i) 
1 — Xx/x'i 



xi g^xi - a?i) 1 

; 2 n / ' i 



7l~ 

12 



xi 



x 



1 L 



(135) 



The calculation of the first term in this expression proceeds identically to the previously 



discussed case. Upon substitution of the second term into the cross section in Eq. (50), 
the delta-functions constraints set X\jx\ = z. Using the relation between hadronic and 
partonic Mandelstam invariants, it can also be shown that they set u = 0. Since the 
Higgs boson transverse momentum can be written in terms of Mandelstam invariants 
as "p\ = tit/Q 2 , the second terms contributes only for zero transverse momentum. The 
final contribution of the X]^ term to the differential cross section is 
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Upon integration over the Mandelstam invariants, we find the following contribution 
to the total cross section: 
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An identical calculation as that outlined above yields the following results for the 
differential and total partonic cross sections contributions from X n : 
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where the differential cross-section term above differs from the corresponding contri- 
butions from X n only by the interchange u ■<->■ t. 

The total next-to-leading order result is obtained by summing the four contributions 
described above: 

d 2 & (l,total) d 2 & (l,a) ^(1,6) ^(l.c) ^(l.d) 
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fractal) = -(!,«) + .(1,6) + .(l,c) + & (l,d)_ (M0) 

We now compare the differential cross section for px > and the total cross section with 



the full QCD expressions given in Eqs. (12) and (13). Summing the four contributions to 
the differential cross section, and keeping only terms which contribute for p? > 0, we obtain 
the following result: 
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We have dropped the plus prescription on the Mandelstam invariants since we are only 
considering the region p? > 0. This expression agrees with the QCD result in Eq. (12) up 
to terms which are finite as u,i — > 0. Such terms can arise from suppressed operators in 
the EFT, and not from emission of collinear or soft gluons. Our result therefore correctly 
reproduces that of QCD in the limit of soft or collinear gluon emission, as expected. 
The total cross section obtained by combining the four contributions above is 
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The terms with 5(1 — z) and ln(l — z) in this result agree with the QCD cross section in 
Eq. ( 13 ). Terms non-singular in the limit z — Y 1 can arise also from emission of hard jets with 
the Higgs, and therefore agreement in the non-singular polynomial in 1 — z is not expected. 
The description of these pieces requires power-suppressed operators in the EFT. Our total 
cross section therefore properly reproduces the expected result from fixed order QCD. 



B. Top down vs bottom up running: impact-parameter space 



Section VI A was devoted to the RG running of the hard Wilson coefficient H(x\X2Q , pi). 
First H(xiX 2 Q 2 , (1q) is obtained from a matching calculation at the hard scale fiQ ~ rrih- 
Next, after computing the anomalous dimension of H(x\x 2 Q 2 , /i), it is RG evolved down 
to the n T ~ p T scale to obtain H(xiX 2 Q 2 , \xq\ \xt) which has the logarithms of m h /p T 
resummed. We refer to this procedure as top-down running. Equivalently, one can compute 
the anomalous dimensions of the iBFs and the inverse soft function and perform an RG 
running from the /it ~ Pt scale, where the iBFs and inverse soft function live, up to the 
hard scale fiQ. We refer to this procedure as bottom-up running. Scale invariance of the 
cross-section implies that the total anomalous dimension of the convolution of the iBFs 
and inverse soft function in the t^^ variables must exactly cancel the anomalous dimension 
of H(xiX 2 Q 2 , /u). In this section, we use this consistency condition [HS1 ES] to obtain the 
anomalous dimension of the iBF which we will compare with what was obtained from the 
fixed order calculation. We will derive the anomalous dimension of the iBF in a third way 



in section VHP as yet another check. 



We define the quantity V as 

V = H( Xl x 2 Q 2 ,fi)X(fi), (143) 
where X(fi) is the convolution of the iBFs and the inverse soft function 



X 



-imax j.— ftnax j- 

S ~ 1{ g^ 6 ^- (144) 



The factorization formula for the differential cross-section at the pt scale in Eq. (41 ) contains 
the quantity T>. Furthermore, all the renormalization scale dependent quantities are con- 
tained in T>. Scale invariance of the cross-section then requires that T> has zero anomalous 
dimension or equivalently that the /^-dependence cancels between H (xix 2 Q 2 , /i) and X(ii). 
We begin by defining the renormalization constants which relate the bare and renormal- 
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ized quantities as 

H(xix 2 Q 2 , fi) 
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where the sunscript b on the objects on the RHS denotes bare quantities and the objects 
on the LHS are the renormalized quantities. Note that the renormalization of the hard 
coefficient of H is multiplicative. On the other hand, the renormalization of the the iBFs 
and the inverse soft function involve a convolution. This implies that the hard coefficient 
has local running while the iBFs and the inverse soft function undergo convolution running. 
The consistency of top-down and bottom-up running implies that the total convolution of 
the iBFs and the inverse soft function undergoes local running that precisely cancels the 
running of H(xiX 2 Q 2 , /i). 

The scale invariance of T> leads to a consistency condition that relates the various renor- 
malization constants as 
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To derive this consistency equation, we first write T> entirely in terms of bare quantities 



and use the first equation in Eq. (145) to replace H^. Next we write T> entirely in terms 



of renormalized quantities and then use the last three equations in Eq. ( 145 ) to replace the 



renormalized iBFs and the iSF in terms of the corresponding bare quantities. Finally by 



equating these two ways of writing T> we arrive at Eq. (146). 



We can make the content of Eq. ( 146 ) more transparent by noting that we can write the 



renormalization constants as a sum of a leading order piece and perturbative corrections as 
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where 



Zs - l( Q> 



148) 



We have used the fact that the dependence on t n /Q and t n /Q of the soft renormalization 
constant is separable [86] as was seen in the explicit calculation in section [V] The consistency 



condition in Eq. (146) can now be written as 
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The leading order results for the soft and hard coefficient renormalization constants from 
sections M and |VI A| respectively are 
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from which we can solve for the one loop renormalization constant of the iBF as 
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which is in agreement with the ultraviolet pole terms for the iBF derived in section |V} 
Thus, we have an explicit demonstration of the equivalence of the top-down and bottom-up 
running at leading order. Eq. (151) implies that the leading order anomalous dimension of 
the iBF is given by 
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In section [Vll D , we provide another independent check by computing the anomalous di- 
mension by taking the derivative with respect to fi of the RHS of Eq. (43) which expresses 
the iBF in terms of a matching coefficient and the standard QCD PDF. 
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C. Top-down vs bottom-up running: momentum space 



The consistency of the top-down and bottom-up running can be checked in yet another 
way by looking at the momentum space factorization. We wrote down the momentum space 



factorization earlier in Eqs. (50) and (53). However, in order to demonstrate the equivalence 



of the top-down and bottom-up running above the pt scale, we need the momentum space 



analog of Eq.(41) which is given by 
d 2 a 



(27T) 



du dt 



i(N t 



dp+dPh / d 2 K / d 2 ki \ d 2 k± \ d 2 kj 5{ki + k£ + ki + kt) 



:2 



5 [u - m 2 h + Qp h ] 5 [t - m 2 h + Qp+] 5 \p+p h - k\ 
J dxi J dx 2 J dt+ J dt^H(x 1 x 2 Q 2 ,^Q]^T) B^(x 1 ,t^,k^,fx T ) B na/3 (x 2 ,t^,k^,fi T ) 



t- t + 
x S(xiQ - pi - x 2 Q -vl - k„, fPr). 



(153) 



The momentum space functions are given in terms of the impact-parameter space functions 
by the Fourier transforms 



B n fi( x i t,k±, fi) 
S(k~,k + ,k ± ) 



d 2 b 



4?r 2 
d 2 b 



- 1 e -ik±-b ± 



47T 2 



—e~ ik± ' b± 



S{k~,k + ,b ± ). 



(154) 



We have made use of Eq. (38) to recast the factorization theorem in terms of purely collinear 



iBFs, defined with a zero-bin subtraction as in Eq. (37), and the soft function. In this mo- 



mentum space version of the factorization formula at the pt scale, the running of the purely 
collinear iBFs and the soft function is determined entirely by scaleless virtual graphs as 
shown in Figs. [4] and |5j The real gluon emission graphs for the momentum space purely 
collinear iBFs and the soft function are finite and do not contribute to the running. This 
is in contrast to the impact parameter space factorization, examined in the previous sec- 
tion, in which the real emission graphs have UV divergences and do contribute to the 
running. The consistency of top-down and bottom-up running should work out in both 
impact-parameter and momentum space. In this section we show how the virtual graphs 
for the purely collinear iBFs and the soft function cancel the running of H(xiX 2 Q 2 , fj) to 
maintain the /^-independence of the differential cross-section. 
For one loop running, we show that 



d 

djj, 



-It? 



-(7g ) + 7g ) V = 7HJ r , 



(155) 
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where J 7 is defined to be the product of the purely collinear iBFs and the soft function 



t- t + 

J= = Bf(x 1 ,t^k^,fiT)B na(3 (x2,t^,k^fi T )S(x 1 Q-p h -^,x 2 Q-pt--^,K^T) 



(156) 



Eq. (155) ensures that Eq. (153) is scale invariant as required and 7^ is the one loop cusp 
anomalous dimension, given in Eq. (114), of the operator O defined in Eq. (92) and 



(i) 1 (i)* 
To +To 



2a s C A 



In 



7T 



A* 



(157) 



The bare and renormalized J 7 are related as 

T h = Z- F F 



T-r, (158) 
and writing the left-hand side in terms of renormalized fields and couplings we get 

ztz 2 



T = B^f Bn a p S + ( 9 — 1) J"' 3 Jna/3 S 



S + 



2{25Z g + SZ A } - bZ T 



B® B na p S 



(159) 



where the quantities B^ , B na /3, S are all written in terms of renormalized fields and cou- 
plings. We suppress the arguments in Eq. (156) for notational simplicity. 
The zeroth order expression for T is 



F (0) ... o(0)a/9 o(0) q(0) 



(160) 



where 



B^(x,t,k ± ^ T ) = -g 2 gfS(t)8(l-x)6^(k A 



S^(k-,k + ,k ± ,ft T ) 



{Nl - 1) 



5(k-)6(k+)6 {2 \k ± ), 



(161) 



and the contribution of the virtual graphs in Figs. [4] and [5] to J 7 is 

jrV(l) = B V(l)a? B ^ S (0) + B WaP gVO) g(0) + g^afi gV^l) 



5V(l)aP r>V(l)aP R (0) Q (0) , 
D n n n0 n na/3 D 

+ B^B%S^\ 



5V(1W r V(1)q/3 
Bn ~ B m 



o(0) or(O) 



(162) 



The one loop virtual graphs contributing to the purely collinear jet functions and the soft 
function are given by 

gV(l)aP _ g(0)a/3 



B 



g(O)a0 



V(l)aP 
nO 

gV(l) = g(0) 



V b (px) + V b (-pi 
Vwipi) + Ho(-Pi) 



V a (pi,P2) + V a (-p 1 ,-p 2 ) 



(163) 
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where expressions for V a {pi,p2) 1 Vb(Pi), and Ho(pi) were given earlier in Eq. (57) and in 
Eq. (103) with off-shell momenta pip to regulate infrared divergences. One can rearrange 
the terms in Eq. (162), at the level of the integrand, using Eq. (163) to get 

_ d(0)(HJd(0) o(0) 



(164) 



x (K(pi,p 2 ) + - Ho(pO] + [H(P2) - Ho(p 2 )] 

+ V a (-pi, -p 2 ) + [Vb(-pi) - Vwi-px)] + [^(-pa) - Ho(-p 2 ) 



The first line in curly brackets is precisely the sum of graphs in Fig. [3] which determine the 
anomalous dimension at one loop for Oin ■ pi,n- P2) and the second line in curly brackets 
is just the sum of conjugate diagrams. From Eqs. (159), (164), (105), and Eq. (110) we get 

a s C A 



1 + 



2tt 



e 2 e \q< 



(165) 



so that 



2a s C A (Q 

7T V jjL- 



(166) 



as expected for the scale invariance of Eq. (153). 



D. Running after matching 

For completeness, we will also verify that the running after matching our results to the 
standard PDFs cancels the running of the hard Wilson coefficient H(Q 2 ,[i 2 ). To do so, we 



begin with the relevant quantity appearing in the factorization formula of Eq. (47): 



J %1 %2 J 



4-max j. j-max f _ ™ ™ 



We have rearranged the x i integrals appearing in Eq. (47) into a form convenient for this 



calculation. The quantity X must obey the renormalization group equation 

dX 2a s C A , Q 2 
fi— = In— X, (168) 

dfJ, 7T jj, 



in order for the cross section to be scale invariant and for the operator O defined in Eq. (92) 
to have the correct anomalous dimension. 

It is simplest to derive the running of X by differentiating directly the beam function 



defined in Eq. (43). This accounts for the running of the product Xx f g / p in Eq. (167). From 
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Eqs. (43), (46), (86), and (91), we can check that the terms appearing in the n-collinear 



beam function that contribute to the \x dependence are 

2/3 



IX J x X \ x l 



5(1 - x[) 



t+ 



Q 2 



xi(l-xi) + 



1 — X-, 



.'L- -i 



+ 



X, 



[1 - x[] 



In 



(169) 



The running of the coupling constant and PDF are given by the standard expressions 



d/j, 

dfg/p(x,fl) 

d\i 



2vr \ 3 3 

dx 



a s f 1 dx fx \ / n a s C A f 1 dx fx \ 



/ , r v 1 — X 

a; (1 — x ) H ; h 



.r 



[1 - x') 



(170) 



We now differentiate the beam function in Eq. (169) with respect to /i, and arrive at the 
result 



dBf 2 a s C 
t 1 j~ = 2 9 ~^9± 



d\i 



2 a s LyA a/3 / 1 dx x I X\ 11,-, , 

" ' -fg/P — ,/X ^(l-XJ 



7T 



xi ^1 



1 



Q 2 



+ tf(f+)tf(l - x'OlB? 

A 4 ' 



(171) 



We note that this result is in agreement with the anomalous dimension found in Eq. (151). 
An identical calculation for the other n-collinear beam function yields 



dB 



a 6 



I'- 



2 Ol s Ca a/3 f l dx 2 f X 2 \ j i, 1 



+ 5(t n )5(l - x 2 )\n% ( . 

(172) 



We next differentiate the finite part of the iSF defined in Eq. (|85|), and find 

ds- 1 



ti- 



dy, 



^i 2 ^v 1 8{ir - o w - tt)^ + ^ 2 m ax - t- f 



Q 2 



j-max f+ 



Q 2 



-/-■max f— 



(173) 



We note that it is possible to manipulate the result of Eq. (85) such that Q rather than Q 



appears in the logarithms and distributions, and we have done so in writing this expression. 
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We are now ready to plug all of this into X and determine its running. The leading order 
expression for X can be easily determined forom the results of Section M to be 



X (0) 



1 dx[ 

—Jg/P 



doc<~) 



X2 ^2 



— fg/P 



X 2 



X 



x g\d-2) 



N 2 



-Q 2 8(1 - x[)8(l - x' 2 )8(tn ax )m 



(174) 



Upon substitution of the derivatives computed above in Eqs. (171), (172), and (173) into 
the derivative of X, we find the following result through next-to-leading order: 



dX 



-2—C A \n— 

7T /i Z 



dx-^ „ 

-Jg/P 



i \ f 1 dx 2 fx 2 

/ —rJg/P[-r>l* 

1 / Jx 2 x 2 \ Uj 2 



T 
X i 



AT 2 - 1 

x g \d - 2) — q 2 5(1 - - siwcwc) 

= -2^C A ln^X(°). 

7T /i 2 



(175) 



The computation is simplified by noting that all plus distribution terms appearing in this 
expression cancel due to the following identity: 



dtt 



8(t; 



-t+ 



8{t n - t r 



o. 



(176) 



such that only the terms in Eqs. ( 171 ), ( 172 ), and (173 ) with an explicit ln(Q 2 / /i 2 ) contribute 
to the running. This verifies that the running of X matches the required form in Eq. (168). 



VIII. CONCLUSION 



We have derived a factorization theorem for the transverse momentum and rapidity distri- 
butions of the Higgs boson at hadronic colliders, in the region m h ^> p T ^> Aq C d, employing 
the methods of effective field theory. The factorization theorem naturally allows for a re- 
summation of the low px region and can be straightforwardly generalized to other processes 
that involve the production of one or more color neutral particles, such as the Drell-Yan 
process. The problem of resummation in the low pt region has been extensively studied 
in the QCD literature [27T429| I3TI 135] and more recently with effective field theory meth- 
ods [4*21 |4"3] . However, the form of our factorization theorem and the method of our analysis 
contains several new and interesting elements that are worth pursuing further. 

The derived factorization theorem takes the form shown in Eq. which is a convolution 
between the perturbative coefficients H(xiX 2 Q 2 , /^q; ^t) and Q^(xi, x[, x 2 , x' 2 ,pr, Y, fix) and 
standard QCD PDFs. This form was obtained by matching onto a sequence of effective field 
theories. First the top quark is integrated out to obtain the standard effective coupling of 
the Higgs to gluons. The hard coefficient H(xiX 2 Q 2 , \iq\ ^t) is obtained by matching the 
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effective Higgs operator onto an effective SCET PT operator, integrating out the hard scale 
ruh in the process, followed by RG evolution between the scales fig ~ rrih and fix ~ Px- The 
RG evolution of the hard coefficient, determined in the effective theory SCET PT , sums up 
logarithms of rrih/px- The perturbative coefficient Q % \ evaluated at the fix ~ Pt scale, takes 
the form in Eq. |2 ) which involves a convolution over the n and n collinear functions 1^ n . g j 
and the inverse Soft Function <S -1 . This form of into two collinear and one soft sectors is 
the result of the soft-collinear decoupling property of the leading order SCET PT Lagrangian. 
The collinear functions X"^ are obtained from an OPE in Aqcd/pt of the impact-parameter 
Beam Functions (iBFs) which appear in the factorization theorem at the px scale as in 
Eq. (41). At leading order in Aqcd/pt the iBFs match onto the standard QCD PDFs with 
Wilson coefficients T^ n . g i as in Eq. (43). The form of the factorization theorem is pictorially 
summarized in Fig. [TJ 

The iBFs are universal objects in that they will appear in the px distributions for any 
gluon initiated process of color neutral particles. The iBFs defined here will have analogs for 
quark-initiated scattering as well, and the generalization to such processes is straightforward. 
Similarly, the iSF is also universal and again has a straightforward analog for quark-initiated 
processes. The iBFs are defined without a soft zero-bin subtraction, so that in addition to 
describing the emissions and virtual effects of collinear gluons they also do the same for soft 
gluons. The presence of two iBFs in the factorization theorem leads to this soft region being 
double-counted. The iSF plays the role of subtracting a soft region so that there is only one 
soft region in the full factorization theorem as required. Each iBF is still defined implicitly 
with an ultrasoft zero-bin subtraction to avoid double counting with ultrasoft modes which 
are unaffected by px constraints. Each iBF and correspondingly the collinear functions 
■^n n-g i are defined in the class of covariant or non-singular gauges. The iBF is not invariant 
under singular gauge transformations where the gauge potential is non-vanishing at infinity 
and must be modified. However, the quantity Q l i in the factorization theorem which is a 
convolution over the two iBFs and an iSF is fully gauge invariant. The factorization of 
into two iBFs and an iSF can be quite useful for the computation of higher order radiative 
corrections at the px scale since the computation of each of the factored objects is much 
simpler. 

The factorization theorem can also be formulated entirely in momentum space with no 
reference to an impact-parameter. This form of the factorization theorem is given in Eq. (fl]) 
with the quantity Q l i given in Eq. (53) entirely in terms of momentum-space objects. In this 
effective field theory formulation one avoids Landau poles associated with impact-parameter 
integrations in the usual CSS formulation. The use of SCET allows power-suppressed cor- 
rections of the form Aqcd/pt and px/rrih to be systematically included. When px ~ Aqcd, 
a model of the non-perturbative region can be included as in the usual approach. 

The structure of our factorization theorem differs also from previous SCET analyses. We 
are led to gluon iBFs that have an index structure, which at higher orders is needed for the 
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separate calculation of effects from each collinear sector. 

We believe our study provides a new framework for the investigation of the low-pr region 
in hadron collider processes. Several interesting differences from previous approaches arise 
naturally in our derivation within SCET that require further investigation. We look forward 
to future research investigating the implications of the factorization theorem derived here. 



Appendix A: Details of deriving factorization in SCET 



PT 



After matching rif — 5 QCD onto SCET PT as in Eq. (27), the differential cross-section in 
the Mandelstam variables u and t is given by 



d 2 a 



1 



du dt 2Q 2 

x 



d 2 Ph , f dn ■ p h dn ■ p h 



(2n)9(n ■ p h + n ■ p h )5(n ■ p h n ■ p h - pf - m 2 h ) 



(2tt) 2 J 2(2tt) 2 

S(u- (p 2 - Ph) 2 )5(t - ( Pl -p h ) 2 ) J2 5^|C(^i,^2)® (hX n X n X s \0{uj 1: uj 2 )\pp) | 2 

initial pols. X 

x (27r) 4 ^ 4 )( Pl + p 2 - P Xn - P Xn - P Xs - p h ), 

(Al) 

where the delta functions in u and t pick out the appropriate final states. The overall factor 
of 1 /4 in square brackets comes from the average over initial proton polarizations. The final 
state hadrons have been broken up into n-collinear, n- collinear, and soft states as 

|X> = \X n ,X n ,X s ), (A2) 

since these are the states with non-zero overlap with the SCET PT operator 0(ui,U2) 



OKwa) = g^Tr\T{S n igB^)^SlS n {gBj^)^Si} 



where 



S n (x^) = Pexp (ig s j dsn- A a s (x» + 

V J oo 



sn") 



(A3) 



(A4) 



with an analogous expression for S n . Since the collinear and soft fields are decoupled from 
each other, we can factorize the matrix element of 0(u>i,u>2) into collinear and soft sectors 
as 



{hX n X n X s \0(uj u uj 2 )\pp) 

(hX n XnX s \g^hTr T{S n {gB^ A _) Un Sl l Sn{gB u nA _) U]n S ] n } \pxp 2 ) 
(X s \Tr{T{S n T A SlS n T B St})\0)(X n \ \ Pl )(X n \ (i^J \p 2 ). 



(A5) 
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The absolute value squared of the matrix element now takes the form 

\(hX n X n X s \C <g) 0\pp)\ 2 
= &®C ® (0|Tr(T{^T D 45 n T c ^})|X s )(X s |Tr(T{^T A 5t^T B 4})|^ 



x 



(Pll (Sg^)^ l^n)(X„| | Pl > 

x (pal (i&i^ |X,)(X,| (fl&J^ |p 2 ), 

= C* (g) C® (0|Tr(f { 1 S fl T z? 4 1 S n T c 4})|X s )(X s |Tr(T{ 1 S n T A 4^T B 4})|0) 

X iV^I (Pl1 ( B WU \ X ^^ ( B ln±a) ui \Pl) 



c 



s 1 

x ]^3i<P2l WW (^ Q ) W2 |P2> 

= C* (g) C(g) (0|Tr(f {^T D 4^nT c 4})|X s )(X s |Tr(T{ 1 S n T c 4^T D 4})|0) 

^n(^i «^ i*»x*»i K^k m 
i 



x 

X 



r (P2| (S^)^ \Xn)(X n \ {B? n±a ) u2 \p 2 ), 

(A6) 



where we have simplified the contraction of color indices in the last two equalities. Next we 
insert the identity operator 

d 4 Pn J d 4 p n I d 4 p s 5^{p n - Px n )t {i \pn ~ Px n )5 {4 \p s - K Xs ), 

(A7) 

and we decompose the momentum components of order m h for the final state collinear 
particles into label and residual parts so that 

p x = Px + K ~x , p x- = Px- + K t, 

(A8) 

where Px , P Xn ~ rrih and K Xn , K x _ <C m h and we write the remaining momentum compo- 
nents as 



Pp 1 - = Kp x , P" 5± = KZ' , P*t = K't . 



(A9) 



Similarly, we write 



Pn — Pn + kn i Pn — Pn + kn > 

n+'- L — h + ' ± rT' 1 - — h - ' 1 - — 



(A10) 
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where again p~,p£ ~ w/i and fc^ <C w^. The delta functions in Eq. (A7) for the large 
collinear momentum components break up into Kronecker delta functions over the hard 
collinear label momenta and residual delta functions which we can write using the integral 
representation to get 

1 = J2 5 p- P <V p+ / d A k n d A k n d%6^(k n - K Xn )6^(k n - K Xn )5^\k s - K Xs ) 



Pn ,pn 



Pn,Pn 



d A x d A y d 4 z 



D i(k n -K Xn )-x i(k n -K X -)-y i{k s -K Xs )-z 



(2tt) 4 (2tt) 4 (2ttY 

(All) 



The exponential factors above involving the momenta of the final state collinear and soft 
states allow us to perform a shift in coordinates and perform a sum over the states X nfitS 
states. Furthermore we break up the Higgs momentum as 

n-p h = n-p h + n-k h , 
n-p h = n-p h + n-k h , 

Pi = (A12) 

where n ■ ph,n ■ ph ~ ran and n ■ kh,n ■ k\, k^ <C m^. The Higgs phase space integrals can 
now be written as a sum over label momenta and integrals over the residual momenta as 

£ J d 2 k h± J d &^6(n .p h + n. p h ) 

x KPtPl + vlK + PhK + KK - k L - m D- 

(A13) 

We also break up the momentum conserving delta function into Kronecker deltas over label 
momenta and delta functions over the residual residual momenta as 

fiWfa +P2 - Pxn - Pxn - Pxs - Ph ) = 5^-5^5^ (K XsX + K Xn± + K Xn± + k h± ) 

x 5(K Xn + K Xn + K Xs + K)5(K Xn + K Xn + k Xs + 

x + K + K + k tMK + K + K + K)- 

(A14) 
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After all of the above manipulations, the SCET PT double differential cross-section can be 
brought into the form 



d 2 a 



(27r) 4 ^ /",,+ ,,_ f ,2k f jii, Mi, Mi, f d 4 y d 4 z 

dudt 16Q 2 (27r)3(iV2 - 1)2 J dk ^ J d ^ J d knd ^ k ' J J^y (2^ (2^ 

x + Ph) S (PhPh + + Ph fc jt + fc /t fc h - 4 2 ± - m 2 h ) 
x <S(u -m^ + Qp fe + QK)5(t -m 2 h + Qp+ + Qk+) 

^ ^ ^ ^ ^ g^fen'^ g^n*?/ gi^s' 2 

initial pols. X n ,Xfi,X s 

x y duidu2 J d{J x dtJ 2 C* {id x , ui' 2 )C(uji, oj-z) 

x W«± a ) wl (o)bi) 
x (pal (s^)^ (v) (^± a ) W2 (o)b 2 ) 

x (0|Tr(f{^T D 4^T c 5t})(^)Tr(T{5 n T c 5t^T D 4})(0)|0) 

x + K + fc+ + fc+)«S(A- + fc" + A- + k~ h ). (A15) 

Next we perform the integral over x + to get the fourier transform of the n-collinear field 
B^ n±/3 (x) with respect to x + as 



/ 



dX + ; ,. , 



e | fe -, (Bf^p) (x + ,x ,x ± ) = (B^p) (k n ,x ,x ± ), 



2(2tt) 

(A16) 



with an analogous equation for the n-collinear field B^^^y). Next we combine the sum 
over labels and integral over residual as usual so that 

J dtuidk~ J duii, J dui2dk£ — > J du>2, (A17) 
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— 4- (') (') 

and absorb the residual momenta k~ and k£ into qj[ and oj 2 respectively. We arrive at 

d 2 a (2ir) f . _ /" , 2l 1 /" „+ ,2,jl /" „- a, 1 ,4, /" dx~d 2 x± dy + d 2 y± d A z 
d^dt = 16Q 2 J dPhdp ~ h J d ^2 7 d *» d *«2 y ^ M *' y 72^l27F (^O 1 
x 6(p+ + Ph )5(p+p h - k 2 ± - m 2 h )5{u -m 2 h + Qp h )5(t -m 2 h + Qp+) 

initial pols. X n ,Xn,X s 

x J du>idu>2 J du' l du' 2 C*{u' 1 ,u' 2 )C{uJi,uj2) 

x (0|Tr(f{5 f - t T D 4^T c ^})(^)Tr(T{5 n T c ^^T D 4})(0)|0) 
x (pi| (Bi ±a )^ (0)|pi> 

x (pal (S&i^ (^±a) W2 (0)|P2> 

x <5 (2) (£; s + fc n± + fc^ + k h± ) 

x 5( W2 - P + - k+ - kt)S(u! -vl-K- K) 

1 

x (N 2 -ir 



(Ms) 



Next we note that the n-collinear matrix element is proportional to 5{ui\ — uj[) which can be 
seen as follows: 

(Pll {BtnLp)^ (*~ *l) (0)|Pl) 

= [< ±/3 (a;-,x ± )5(pt _ ^ |X n )(X n |5(P - c^X^O)^) 

= E <^!| [^^(^-^-L)^^! - ^i)] |^><^| [^(^ - ^ 1 ) J B^ ct (0)] |^>, 
= (Pll [^(s" - ^)] - Wl )<Xa(0)] bl), 

(A19) 

where we have inserted a complete set of states and used momentum conservation and 
then performed the sum overs states again. Analogously, the n-collinear matrix element is 
proportional to 5(u 2 — oj 2 ) so that one can perform the integral over u[ 2 to get 
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d 2 a 2tt 
dudi ~ 16Q 2 (A 2 - l) 2 

dx~d 2 x± dy + d 2 y± d A z f db + db~ d 2 b 



d Pi d Ph J d 2 h ± ^ J dk+d 2 k^ J dk n d 2 k±d 4 k s 



x 



(2tt) 3 (2tt) 3 (2tt) 4 J 2(2vr) 2 (2tt) 2 
x 0(p£ + P h )KptPh ~ 4x - - m 2 h + Q Ph )5{t -m 2 h + Qp+) 

x e |^(a; _ -fe _ ) e -^nX-(^±-bx) e %k^(y+-b+) e -ik fl± -(y ± -b ± ) 
X e l fc ^ (^~- 6 ~) e f fc ^(^ + - b+ ) e -*^ s ±-(2±-bx) e ^; i x-6x 

(A20) 

where we have defined the jet and soft functions 

initial pols. 

■tfW+yj.) = J2 (P2\[9Btn±,(y\y±)S(V-u J2 )gBi ±a (0)]\p 2 ) 

initial pols. 

S(z) = (0|Tr(f{^T D 45 n T c ^})(^)Tr(T{5 n T c 5t^T D 4})(0)|0). 

(A21) 

After performing the integrals over the residual momenta k+, k^, k^, k^, k^ and the x, y, z 
coordinates we arrive at the simpler form 

2 " f dpUpT I d 2 k h± [ du.du, ! ** e ^-^e^^^ e fa* ± 



dudt ~ 8Q 2 (A 2 -1) 2 J Fh Fh J h± J J 4 (2vr)4 

x d(pt +Ph)$(ptPh - 4 2 ± " m 2 h )S(u -m 2 h + Qp h )5(t -m 2 h + Qp+) 
x \C(u x , u 2 )\ 2 Jfiu!, b-,b ± )jffa, b + , b ± )S(b + , b-, b ± ), 

(A22) 



which appears in Eq. (34). 



Appendix B: Equivalence of zero-bin and soft subtractions 



In this section we demonstrate the validity of Eq. (38), which we write here again for 
convenience 

E = J du 1 du 2 \C(uJi,uj2,[j 1 )\ 2 J dk+dkaB% p (u u k+,b ± ,fj.) B na(3 (uj 2 ,k^,b ± , n) 
x S(u 1 -p h - k n ,u 2 -p+ - k+,b ± ,(j) 

= J dwidw 2 \C(ui,u 2 , /i)| 2 J dk+dk^B^(uJ 1 ,k+,b ± ,n) B na/3 (uj 2 ,k^,b±, p) 

X S^Wl ~Ph ~ k n^2 -Pt - &_!_,//). 

(Bl) 
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The statement of Eq. (Bl) is that the factorization theorem of SCET PT in terms of the 
purely collinear iBFs B^%, defined with a zero-bin subtraction, and the soft function S can 
be equivalently formulated in terms of the iBFs B^ n , defined without a zero-bin subtraction, 
and an iSF S^ 1 . At zeroth order the purely collinear iBFs B^^ 13 and the soft function are 
given by 

B^(tu u k+,b ± ) = -g 2 gf6(k+)6(u 1 -n- Pl ), 



B n (^2, K n > b -i 

5 (0) (wi -pi -k^,u 2 -p+- k+,b A 



-9 9± P 5(k n )5(w 2 ~n-p 2 ), 
(N 2 -l) 

'-5 (ux - Ph - k n )8(u 2 - vl - k+ _ 



(B2) 



and Eq. (Bl) is trivially satisfied since 

B^(u h k+,h 
B^(oo 2 ,k^,b A 
S (0) (uj 1 -p h - k n ,u 2 -p+ - k+,b A 



B^(u u k + n ,b ± ) 

S^-^-pl-kz^-pt-k+M 



(B3) 



We now demonstrate the validity of Eq. (Bl ) at the next order in perturbation theory. The 
the purely collinear iBFs and the soft function up to the first order in perturbation theory 
take the form 



Bf{u x XM = B^(u l7 k+,b ± ) + [BWMfaXM 
+ ■■■ 

Bf(u 2 ,k^,b ± ) = B^(u; 2 ,kT,b ± )+[B^(co 2 ,kT,b ± ) 
+ ■■■ 

S&i-Pl -k^,u 2 -p£ -k+,b±) = 5 (0) (wi-Pfc -k^,u 2 -p+ -k+,b ± ) 

+ S [1 \ui-p h -K,uj 2 -pl-k+,b 1 _) + 



where the dots denote terms that are higher order in perturbation theory and B 



(l);a/3 
nOMO 



de- 



note the zero-bin subtraction terms. There are two types types of terms at first order in 



perturbation theory for the convolution of the left-hand side of Eq. (Bl) so that we can 
write 



E = E<® + + • • • 



= E 



(i) 



E 



(i) 

B i 



(B5) 
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where is the zeroth order term, is the first order term, and E A and E B denote 
the two types of terms at first order in perturbation theory and are defined as 



E$> = j dwidw 2 |C(wi,W2)| 2 J dk+dk: 



'flf^K k+, b ± )B^(co 2 , k^, b ± )S®{m -p^ + k^,u 2 -pt + k+, b ± ) 



(B6) 



and 



£^ ee J duidu 2 \C(uJi,uj 2 )\ 2 J dk^dk~ fi 



x 



■si ^KA^,6j.)S^( Waj A^ l 6j L )5«(wi -Pfc - K,u 2 -p+ - k+,b ± ) 
Bf^^M^iU^' fc«.6±)5 (0) (wi -Pfc - k^,u 2 -pt - KM 



(B7) 



At the first order in perturbation theory, Eq. (Bl ) is equivalent to demonstrating that E 
in Eq. (B7) can be rewritten as 

dcuiduj^C^,^)] 2 I dk+dk~ ft 



(i) 



E 



(i) 



(B8) 

which is the statement that each of the zero-bin subtraction terms in the last two lines of 



Eq. (B7) is equivalent to a soft subtraction. Using the zeroth order expressions for the iBFs 
and the soft function from Eq. (B2) in Eq. (B7) we can write Eg as 



E 



(i) 

B 



X 
+ 
+ 



duidu 2 \C(ui,u 2 )\ 2 

2g i S(ui - n ■ pi)5(u 2 - n ■ p 2 )S {1) (u! ~Ph,u 2 -p%,b±) 
9 2 (N 2 - 1) 



4 

9 2 (N 2 C - 1) 



5{u 2 - n ■ p 2 )5{u x - p h )g^B^. afi (u 1 , u 2 - p+ b A 
5{ux - n • j?i)<y(w2 - pt)gfB^ Ju 2 , ui - fej 



(B9) 
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The iBFs, the zero-bin subtraction terms, and the soft function at first order in perturbation 
theory receive contributions from virtual (V) and real (R) emission graphs so that 



B 
B 



(l);a/3 
n,n 

(l);a/3 
nOMO 



B 



+ B 



S (D = S V(1) +S R(1) 



and correspondingly we write 



E 



(i) 

B 



E 



V(l) 
B 



+ E 



H(i) 

B ' 



(BIO) 



(Bll) 



where E^ 1 ^ and E^^ denote the contributions from virtual and real graphs respectively to 
E^ and are given by 

E V z R{1) = J du: l dw 2 \C{uj l ,uj 2 )\ 2 

x 2g i 5(u 1 - n-p 1 )5(u 2 -n-p 2 )S VMl \uj 1 - j£, u 2 - p£, b±) 
g 2 (N 2 — l) a p v,R(i), + , \ 

+ 4 <H W 2 - n • p 2 )S(u 1 - V h )9± B nO;ap ( Ul i U2 ~ Ph> 6 -U 

g 2 (N 2 - 1) _ + a/3 v,R(i)( — j \ 

+ ^ o(^-n-pi)d(w2-p^/S^.^(w 2 ,wi-p h ,6 ± ) . 

(B12) 



We demonstrate the validity of Eq. (Bl ) at first order in perturbation theory, or equivalently 



Eq. (B8), separately for the virtual and real contributions. We first look at the virtual 



contributions for which the soft function and the zero-bin subtraction terms of the purely 
collinear iBFs take the form 



B 



V(l)aP t 
nO, 



V, cu 2 -p+,b ± ) = B^^, u 2 -p + h ^){-2ig 2 C A )I s 
5 Io 1)Q/3 ( w 2, u 1 -p h ,b 1 _) = J^ 0)af3 (uj 2 , u 1 -p^,b ± )(-2ig 2 C A )I s 
S v U(u 1 -pJ t ,u; 2 - P +,b 1 _) = - Ph ,u 2 - pt,b ± )(-2 ig 2 C A )I s 



(B13) 



where I s is the scaleless integral 



d d 



{27r) d {£ 2 +iO)(n-£- iO) {n-i + iO)' 



(B14) 
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Using Eqs. ([B13J), (jB14j), (|B2j), and (|B3j) in Eq. (|B12j) for we find that 

E V B {1) = -2g A J duj l( iuj 2 \C(uj u uj 2 )\ 2 

x 5{ui - n-pi)8(u 2 - n ■ p 2 )S v{1 \uj 1 -p h ,u 2 -p%,b±) 
= — J duiduj 2 \C(uji,uj 2 )\ 2 J dk+dk^ 



(B15) 



which is in agreement with Eq. (B8) for the virtual graph contributions at first order in 
perturbation theory. 

Next we look at the real contribution E B for which we need expressions for the single 
gluon emission contribution to the iBFs and the soft function which are given by 

J4i 



s^-p-^-ptM = g 2 c A {Nl-i) J j^e(k°)6(k> 



x 8(ui —p h - k )5(u 2 -pl~k + ) 



k+k~ 



B 



i?(l);a/3 
nO 1 



uJi,u 2 -pi,b±) = -gjAg C A 



d 4 k 



6(k°)5{k 2 )5{n-p 1 -uj 1 -k-) 



x 5{u 2 -p+ - k + ) 



B 



no 



k+k~ 

'4/ 



Kwi-Pfc.fci) = -9l^9 4 C A J A^(2n)6(k )5(k 2 )8(n.p 2 -u 2 -k 



x S(ui -p h - k ) 



k+k~ 



(B16) 



Using the expressions in Eq. (B16) in Eq. (B12) for E B and using the freedom to perform 
residual shifts in the label momenta U)i j2 , but ignoring these shifts in the Wilson coefficient 
C(uj\,uj 2 ) as they are power suppressed, we find 

e b W = ~ 2 9 4 J duidu 2 \C{ui,uj 2 )\ 2 5{uji -n ■ pi)5(u 2 -n-p 2 )S {l \ui - pf^,u 2 - p^,b±) 

duj 1 dtu 2 \C(uji,uj 2 )\ 2 J dk^dk^ 



x 5J>) ^(u^, b^mj^KMS^i^ -pX,u 2 -p\ - kt,b ± ), 



(B17) 



which is in agreement with Eq. (B8) for the real graph contributions. Thus, together with 
Eqs. flB"l5] ), ( |B17fr , p8| ), ( |BTo| ), and pll) we have explicitly shown the validity of Eq. flBl| 
up to first order in perturbation theory. 
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Appendix C: Consistency of top-down and bottom-up running above 



In this section we give details of the derivation of the consistency relation in Eq. (146). 
After performing a shift in the integration variables, the scale invariant quantity T> 

V = J (lt n J dtnH(x 1 X 2 Q 2 B^(x 1 ,t n ,b±,IJ 1 ) Bnaf}(x2,tn,b±, /j) 
xraax j._ xmax x 



Q 



Q 



(Ci) 



can be written as 



dt n J dt n H( Xl x 2 Q\ n) Bf( Xl ,t n + t™ ax , b±, fi) B naP {x 2 , t n + Q™, b ±} //) 
xS (--,--,6^). 



(C2) 



We now write the above expression in two equivalent ways. First we write T> entirely in 



terms of bare quantities and then use the first relation in Eq. ( 145 ) to write the bare hard 



coefficient Hf, in terms of the renormalized hard coefficient H to get 



V 



K J K J K J < h{x iX2 q\ p)z H 5(t" n -t' n - tz x Wk -t'n- 1\ 



(C3) 



Next we write D entirely in terms of renormalized quantities and then use the last three 



relations in Eq. (145) to write the bare iBFs and the bare iSF in terms of renormalized 
quantities to get 

V = ^ J < J < J dt' n J dt'J J ds n J ds n J dt n J dt n H(x 1 x 2 Q 2 ,fi) 

—t- 4-f' — t 4- t' 
X Zj n {S n — Z n , \l)/j n \S n — l n , /J,)Zj s _ 1 { — , — ,H) 



x B^x.X^ri Bn, b -M^ 



(C4) 



Since the expressions on the RHS of Eqs. (C3) and (C4) must be equal, by comparing the 



quantity in curly brackets in Eq. (C4) with the integrand of Eq. (C3) one can derive the 
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relation 



Z H (p)5{s n -t n )5(sn-tn) = ^ J dt' n J d4Z n 1 (^-s n ,//)Z n 1 (4-s fi ,//) 

-1 / ~b ^fi ~b ^r< 



as quoted in Eq. (146). 



x Z-\( 



Q 



Q 



(C5) 
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